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SUMMARY 


A statistical model of drop-impact erosion rate has been developed 
in terms of occurrences at the site of a typical cell. The model is re- 
stricted to brittle materials whose strength and/or energy-absorbing 
capability is small in comparison with the energy delivered by each drop 
impact. Consideration of the characteristic mode of failure of solids 
under single-drop impact (formation of a circular crack) leads to the 
result that, if each impact produces fracture of the solid struck, it 
is impossible to remove one fragment from each cell in an array of cells 
without removing a layer of material that is many fragments thick. 

The effect of changing angle of attack as surface roughness increases 
with test time is reflected in the number of impacts required to remove 
an eroded fragment; the required number of impacts includes non-crack- 
forming impacts which result only in work-hardening if the solid under 
test is a metal. 

An expression for the rate of fragment ejection has been developed 
by summing the binomial probability for x hits in n trials between 
limits imposed by a counting rule. The counting rule is a function of 
relative impact velocity, drop mass, strength of the solid, angle of 
attack of the impinging drops, work-hardening capacity of a metal, and 
accumulated subsurface damage. The volume of a typical fragment has been 
found to vary as the third power of the impact velocity if the eroded 
fragment is a tetrahedron and as the fourth power of the impact velocity 
if the eroded fragment is a prism. This velocity dependence is reflected 
in the rate of volume loss. 

A partial test of the rate equations has been performed using drop- 
impact volume-loss data for five selected metals. The agreement between 
theoretical and experimental volume loss is acceptable when the theo- 
retical loss is obtained with eroded fragments of equal (average) size. 

The agreement is substantially improved when fragments of increasing 
volume are considered. 

Correlation has been found between the 0.2 percent offset yield 
strength and the number of impacts required to remove an eroded fragment. 

xi 


Correlation has also been found between the ultimate tensile strength and 
both the volume of an eroded layer and the percentage change in layer 
volume. 

The statistical model is not completely developed since two quantities 
that are needed to calculate the rate of erosion with use of the equation 
that is given must be assessed from experimental data. 

For impact velocities of 1000 ft/sec or less, materials having 
strengths comparable to that of aged Udimet 700 alloy are beyond the 
limits of applicability of the statistical model in its present state 
of development. 
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A MODEL FOR MULTIPLE-DROP- IMPACT EROSION 
OF BRITTLE SOLIDS 


I. INTRODUCTION 

When liquid drops or jets impinge against the planar surface of a soiled, 

high pressures and high fluid-flow velocities are produced. If impingement 

is continuous, even the strongest materials eventually fail. Cavitation 

(l) a 

damage is similar to drop-impact damage. Naude and Ellis observed that 

a high-speed jet forms during the collapse of a cavitation bubble. This 
observation has provided a basis for understanding the similarities that 
exist . 

Eventual development of a theoretical expression for the rate of drop- 
impact and cavitation erosion is important for several reasons. From an 
applied standpoint, such an expression would make it possible to assess 
the probable weight loss from an eroding part over an arbitrary time in 
service. It would also provide a basis for rating the resistance of ma- 
terials to this form of attack. In addition, the very process of developing 
and testing such an expression should produce valuable insights into the 
mechanism of the erosion process. 

The average rates at which cavitation and drop-impact erosion progress 

have been assessed from experimental weight-loss data for a number of metals 

( 2 ) 

by several groups of investigators. Thiruvengadam and co-workers plotted 
an average rate against total test time. They found curves with the follow- 
ing features: a period during which essentially no weight loss is observed 
(incubation period) , a period characterized by a rapid increase in erosion 
rate (accumulation period), a maximum rate, a decrease in rate (attenuation 
period), and a steady-state rate. Some variations of these features have 
been observed by other investigators. Hammitt, Robinson, Siebert, and 

( 3 ) 

Aydinmakine found two and even three maxima in some but not all of 
their plots. 


lumbers in parentheses refer to literature references at the end of 
this report. 
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( 4 ) 

Heymann discussed time dependence of erosion rate in terms of life- 
times of layers of cells which were described by statistical distribution 
functions; numerical results were obtained by assuming that no damage is 
done to underlayer material until it has emerged to the surface. Mok, 
on the basis of the fact that damage to underlayer material does occur if 
stresses within this material exceed a critical value, produced a treatment 

of erosion rate based on cumulative fatigue damage; his treatment does not 

( 2 ) 

give the initial or final features of experimentally determined curves. 


Hammitt, Huang, Kling, Mitchell and Solomon*^ developed a treatment 
of erosion rate based on the energy-flux concent of Hoff, Langbein, and 
Rieger^^ ; numerical results were obtained by excluding materials of high 
erosion resistance such as the tool steels and highly resistant alloys 
including the family of Stellites . 

.O') 

It was suggested' that a treatment of erosion rate based on lifetimes 

of individual cells might prove to be informative because it would lead to 

consideration of successive states or conditions that exist as erosion attack 

progresses. As a first approach to such a model of erosion rate, a simple 

/O') 

but physically unreal case was considered. ' 


Lifetimes of the members of an assembly of cells can be treated sta- 
tistically by considering the lifetime of a cell in the original surface 
layer and the lifetimes of successive underlayer cells below it through 
the thickness of the test specimen. 1 In the first development ot this 
cellular model (see Section II) , the erosion process is idealized by assum- 
ing that all of the eroded fragments have equal (average) volume. The 
attractive feature of this simplifying assumption is that it makes it 
possible to calculate rate of erosion, to calculate accumulated volume 

loss from this rate, and to compare the calculated volume loss with that 

( 9 ) 

obtained experimentally. However, it has been reported by Gould that 
the size of eroded fragments increases with lapse of test time. A second 
development of the cellular model (see Section IV) takes the feature of 
increasing eroded fragment size into account. 

^1 am indebted to Dr. Wayne Nelson, Information Sciences Laboratory, 

General Electric Company, Schenectady, New York, for the suggestion that 
the cellular model can be developed in terms of a typical cell. 
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II . ST ATISTICAL MODEL BASED ON ERODED FRAGMENTS OF EQUAL VOLUME 

Development of an equation for volume rate of erosion, in E of this 
section, is preceded by a consideration of major factors that affect this 
rate, in A through D of this section. 

A. Strength of the Solid Under Test 

The characteristic appearance of surface roughness that develops as 
drop-impact erosion progresses is different for materials of different 
strength and different ability to accept energy elastically. The extreme 
cases are: (A) materials whose ability to accept energy without fracture 
is small in comparison with the impact energy delivered by each liquid- 
drop blow, and (B) materials whose ability to accept energy without frac- 
ture is large in comparison with the impact energy delivered by each liquid- 
drop blow. For any given material, the classification is relative; it may 
be either in Class (A) or in Class (B) depending on the impact velocity and 
drop mass used in testing it. 

The surface of a brittle material of Class (A) is uniformly eroded. 

For a material in this class, it appears that long-term drop-impact erosion 
may involve movement of an eroded surface layer through the thickness of 
the test specimen (see Figure 1,A). A Class (B) material will start to fail 
at weak spots; these spots may be grains with an unfavorable orientation or 
points where several grains meet as was suggested by Von Schwarz and 
Mantel . When eroded fragments are ejected from the weak spots, residual 

crack ends remain. Residual crack ends are able to grow in length as more 
drop impacts occur with the result that eroded fragments will continue to 
be ejected from the weak spots. Pits eventually form at the sites of the 
weak spots and deepen until they pierce the test plate (see Figure 1,B). 

The model of erosion rate in drop impact and cavitation developed in 
this report is directed to materials of Class (A) . This means that each 
impact does some damage to the solid struck (work-hardening or fracture) . 

B. Geometry of Fracture 

Baker, Jolliffe, and Pearson^^ have observed that, for corresponding 
positions on curves of mass loss against mass of impinging water, surface 
coarseness increases with drop size and distance between adjacent erosion 
peaks is proportional to and of the same order as the drop diameter. It 
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A. Brittle Material of Low Strength 



B. Brittle Material of High Strength 


Figure 1. Progress of Erosion on Two Brittle Materials of Widely 
Different Strengths. 


4 




( 12 ) 

has also been observed that the fracture pattern produced by drop 
impact is a circle or polygon of cracks depending on whether the solid 
material is isotropic or anisotropic; distance across this fracture pat- 
tern is roughly equivalent to the radius of the drop that impinged. Fig- 
ure 2 shows such a circle of cracks in hot-pressed alumina. 

The geometry of fracture appears to be dictated by the way in which 
an impinging drop produces stresses in tension and shear in a solid. When 
a liquid drop impinges normally against a planar solid, the impact pressure 
exists over a circular area. By analogy with the impact of a solid sphere, 
radial tensile stresses around the compressed area must be imposed and 
released during the impact. If the tensile strength of a brittle solid is 
exceeded by these radial tensile stresses , a single drop impact will pro- 
duce not just one crack but all the cracks that comprise the circular (or 
polygonal) fracture pattern. In order to develop a cellular model of drop- 
impact erosion rate, it is necessary to assess the number of cells that will 
receive cracks during a single impact because the circular fraction pattern 
passes through them. 

1. Number of Cracks per Impact Produced in a-Cells 

Let drops, all of which have the same radius, r, impinge at normal 
incidence over an area A of a test specimen and let the area A consist of 
an array of cells all of which have the same area, a. Designating the 
number of cells in area A with the symbol q, 

n = A/a. (1) 

Let the array of a-cells be represented by a matrix of closest packed 
circles; let the cell diameter be d. See Figure 3. Let the fracture 
pattern be a circle having radius r/2 and let the circle of fracture be 
inscribed on the matrix of a-cells. A certain number of a-cells is cut 
by the inscribed circle of fracture; this number is designated by the 
symbol n 1 • 

Three circles of fracture with radii Id, 2d, and 3d, respectively, 
have been inscribed on the matrix of a-cells shown in Figure 3. At the 
bottom of Figure 3 the count of a-cells cut by the inscribed circles of 
fracture is tabulated along with corresponding numbers obtained by 
dividing the circumference of the circle of fracture by the cell 
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1 mm 



(a) Projectile plate: hot-pressed alumina; 
target: 2-mm mercury drop; impact 
velocity: 1403 ft per sec. 



(b) Schematic Drawing of Subsurface 
Damage 


Figure 2. Damage Mark Produced on Leading Face of Specimen of a Brittle 

Material as a Result of Collision with a Liquid Drop. View (a) 
is Reproduced from Reference [12]. 
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Number of o-Cells Cut 


Radius of Circular Crack, r/2 Actual Count wr/d 

d/2 0 3.1 

d 6 6.3 

2d 12 12.6 

3d 18 18.8 


Figure 3. Estimate of Number of o-Cells Cut by the Circular Fracture Pattern. 
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diameter, d. It can be seen that the extent of agreement that exists is 
adequate to justify the relation 


n' = rr/d (2) 

as an estimate of the number of cells that develop cracks when a single 
drop impact occurs. 

In the cellular treatment of erosion rate that follows, an eroded 
fragment is identified with _a single a - cell and the diameter . _d, of _a 
g -cell is taken to be the same as the distance across an eroded fragment . 

2. Cracks Produced in a Typical Cell 

In the special case illustrated in Figure 4, the center of impact 
occurred at the center of cr-cell W*. This is also the center of the 
circle of fracture which contains V cr-cells. Let cr-cell W be regarded 
as a typical cell. If the centers of three drops that impinge against 
the circle of fracture shown in Figure 4 fall in three of the cr-cells A* 
through V , the three circles of cracks formed as a consequence of these 
impacts will cut the typical cell W . In Figure 4 the cr-cells struck 
have been taken to be A*, I*, and P*; the center of impact in each of 
these cr-cells is indicated with a solid circle. Each of the circles of 
fracture generated by these impacts will also cut through a total of 
V cr-cells; for each of these three impacts, the typical cell W will 
be only one of the V* cr-cells cut. 

The probability of hits in the typical cell W* is statistically the 
same as the probability of hits in any of the cr-cells in the array in- 
cluding cr-cell W* and cr-cells A* through V*. Physically, the probability 

5 ^* 

of hits in cr-cells A through V determines the probability of cracks in 
cr-cell W whereas the probability of hits in cr-cell W determines the 
probability of compression in this cell. Because these probabilities 
are identical, the probability of hits in cr-cell W determines not only 
the probability of compression in cr-cell W but also the probability of 
crack formation in cr-cell W . In short, the probability of crack for- 
mation in cr-cell W is a function of hits in cr-cell W but the cracks 
that form in cr-cell W* are not caused by hits that occur in cr-cell W*. 
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Figure 4. Schematic View of Macrocell Consisting of V* o-Cells 
Surrounding a-Cell W*. (n' = V*) 
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3. Fragments Ejected from a Typical Cell Site 

Reasoning based on the preceding considerations leads to the number 
of cracks produced in a typical cell site when each of the o-cells in an 
array of such cells has received one or any other arbitrary number of 
impacts. Let the condition on the rain of drops be that it is a random 
rain in which the probability of impacts is uniform. For the purpose of 
the argument developed in this section only, let the further restriction 
be imposed that the distribution of impacts is uniform not in the usual 
sense that each cell has an equal probability of being hit but that each 
cell is actually hit the same number of times. 

The number of a-cells in the circle of fracture is V. If each 
drop impact produces n' cracks, and if the specified type of rain pro- 
gresses until one impact has occurred in each of the n cells in the 
impingement area A, then the total number of cracks formed is nn T and 
the number of cracks formed per a-cell is n'. Similarly, if the rain 
progresses until two impacts have occurred in each of the n cells, then 
a total of 2nn' cracks, or 2n* cracks per cell, have been formed. And 
if the rain progresses until three impacts have occurred in each of the 
n cells, then a total of 3rm' cracks, or 3n' cracks per cell, have been 
formed. 

The number of cracks formed is an integral multiple of the product 
nn'. A rigorous proof of this statement has not yet been developed but 
the following intuitive argument indicates that it is correct. Because 
it is impossible to deliver a fractional part of an impact in each of 
the n cells, the only way that a fractional part of n' cracks could be 
formed per cell is that some cells would receive more impacts than others. 
This would require a nonuniform distribution of impacts which would 
violate the restrictive condition that the random rain has a uniform 
distribution of impacts. 

Let it be required provisionally that N Q cracks must form in a. 
a -cell before an eroded fragment will be ej ected and let the rain pro- 
gress until each of the n cells in the impingement area A has received 
N Q impacts. Let one of the n cells be selected as a typical cell. This 
cell is surrounded by n’ cells in the circle of cracking. If each cell 
in the array has received N 0 impacts, then each of the n' cells in the 
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circle of cracking around the typical cell has received N Q impacts and 
each of these impacts has produced one crack in the typical cell. If 
all of the cells in the array have received N 0 impacts, the typical cell 
has received N q impacts and, in addition, the typical cell has received 
N 0 p* cracks. Because it was specified provisionally that the formation 
of N Q cracks in a cell results in the ejection of a fragment, the site 
of the typical cell must eject p' fragments. Because there are p cells 
in the impingement area A, the total number of fragments ejected from 
all of the cell sites will be nn * . 

Because of the geometry of liquid-drop-impact fracture, it is not 
possible for each cell in an array of cells to lose one and only one 
fragment. It is possible for each cell to lose p* or some integral 
multiple of p' fragments. Although the argument that excludes all but 
integral multiples of the product pp' is so far supported only by an 
intuitive proof, it is confirmed jt posteriori by experimental evidence 
(See Section III .1) . Because of the geometry of liquid-drop- impact 
fracture , it is not possible to remove one fragment from each of the 
g -cells in an array of g -cells without simultaneously removing _a layer 
of material from the test specimen ; the layer removed consists of pp' 
fragments. 

C. Angle of Attack 

For impacts that occur at angles smaller than 90 degrees, the amount 
of energy that is delivered to the material of the test specimen by each 
drop blow is reduced. Erosion investigators have compensated for this 
effect by considering only the normal component of the impact velocity. 

The normal component of the impact velocity, V is V sin 0, where 0 is 
the angle of attack measured from the surface of the test specimen. 

In terms of the normal component of the impact velocity, the energy 

2 

delivered by impinging drops is 1 M (V sin 0) where M is the mass 

2 

of a drop. If impingement occurs at normal incidence against the original 

2 

surface of the test specimen, the energy delivered per impact is JL MV 

2 

because the sine of 90 degrees is unity. Let us identify this energy 
per impact with the formation of cracks of a given length in a particular 
solid. 
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After the first layer of nh* fragments has been removed from the 
original, polished surface of a test specimen, a rough surface of under- 
layer material is exposed; the drops now impinge against this rough 
surface. Without the guidance of specific information, let it be assumed 
that on the average the angle of attack has been reduced to 45 degrees. 
Then the energy delivered per impact will be halved because the sine of 
45 degrees is VF /2. If the energy delivered per impact is halved, 
and if the fractional part of the impact energy invested in crack for- 
mation remains the same, the surface area of the cracks that form must 
also be halved. 

This suggests that twice as many crack-producing impacts may be 
required to circumscribe a fragment when the angle of attack is 45 degrees 
than were required when the angle of attack was 90 degrees. In Section 
II. 3, the number of impacts required to circumscribe a fragment with 
impingement occurring at normal incidence was provisionally taken to be 
N q . If the number of impacts required to circumscribe a fragment when 
the angle of attack is 45 degrees is N-^, then, on the basis of the pre- 
ceding reasoning, N-^ ^ 2N Q . 

After the loss of the second layer of fragments is an accomplished 
fact, the angle that the exposed underlayer material presents to the 
impinging drops may or may not change again for a given material. If 
angle of attack is the only factor affecting the number of crack-producing 
impacts required to circumscribe a fragment at this stage of the erosion 
process, recourse to experimental data (See Section III) suggests that 
the angle of attack must be reduced to about 30 degrees. At this angle 
of attack, the energy delivered per impact will only be one fourth as 
much as it was when the drops impinged at normal incidence because the 
sine of 30 degrees is 1/2. The surface area of the cracks produced at 
this angle of attack will only be one fourth as large as it was when 
impingement occurred at normal incidence and, if the fractional part of 
the impact energy invested in crack formation remains the same, four 
times as many impacts should be required to circumscribe the fragments of 
the third layer of material to be removed from the test specimen as were 
required in the case of the first layer. If the number of impacts re- 
quired to circumscribe a fragment when the angle of attack is 30 degrees 
is N 2 , then, on the basis of the preceding reasoning, N 2 = 4N C . 
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The number of impacts required to circumscribe a fragment after a 
given number of impacts has already been sustained may be represented by 
the sequence 


N i= V N l’ N 2 


(3) 


which is hereafter referred to as the counting rule . The evaluations of 
N Q , N^, and ^ that were made in the preceding paragraphs are based on 
the premise that angle of attack governs the magnitude of successive 
counting-rule numbers and that angle of attack is progressively reduced. 
However, angle of attack is only one of a number of variables which 
determine the magnitudes of the counting-rule numbers. 

D. Work-Hardening of Metals 

Only a brittle metal, whose fracture strength is exceeded when it 
sustains an impact at an arbitrary velocity will develop a crack as the 
result of this impact. Most metals are ductile and must be work-hardened 
to the point of embrittlement before a crack will form. A question that 
needs to be answered is : How many superposed liquid drop impacts are 

required to work-harden a given metal to the point of embrittlement? 

The smallest possible number would be just one and, although it is not 
conclusive, the following empirical evidence suggests that one impact at 
a high velocity may be adequate for some metals. 

Knoop microhardness numbers taken through the thickness of an alclad 

aluminum alloy specimen subjected to waterdrop impacts at the Cornell 

Aeronautical Laboratory, Buffalo, New York, shifted from a maximum of 

74 at the eroded surface to an asymptotic value of 55, which was the 

(13) 

characteristic hardness of the metal. ' Although the test velocity 
was not specified, it was very likely 805 km/hr (500 mi/hr) because this 
velocity was commonly used in performing tests with use of the rotor and 
artificial rain facility at this laboratory at the time that the speci- 
men was tested. 

The ratio of maximum hardness produced by the rain impacts to 
asymptotic hardness of the metal is 74/55 or 1.34. Two conclusions can 
be drawn from this information: first, multiple waterdrop impacts at 

the specified velocity were able to produce a 34 percent increase in the 
hardness of the specimen, and, secondly, because the specimen was severely 
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eroded, a 34 percent increase in hardness can be associated with a 
sufficient embrittlement of the alclad aluminum alloy specimen material 
to permit crack formation in it. 

Knoop microhardness numbers were also taken around a single crater 
produced in 1100-0 aluminum by the impact of a steel sphere at a velocity 
of 626 km/hr (389 mi/hr) . The average ratio of maximum hardness to 

asymptotic hardness for measurements taken below the crater and from the 
edge of the crater near the surface was found to be 1.54. This means 
that a 54 percent increase in hardness was produced by a single steel- 
sphere impact. This is in excess of the 34 percent increase in hardness 
required for sufficient embrittlement of the metal to permit drop-impact 
erosion of an alclad aluminum alloy specimen to progress. 

Craters produced in 1100-0 aluminum by single waterdrop impacts at 
velocities above 1000 ft/sec were also studied but, regrettably, 

Knoop microhardness numbers in the metal around these craters were not 
taken and the tested specimens are no longer available. Consequently, 
it is not known if a single waterdrop impact at these high velocities 
is able to produce work-hardening to the point of embrittlement of metals 
such as 1100-0 aluminum. 

On the basis of the preceding evidence, it can be assumed that, for 
the high velocities being used in applications where drop-impact erosion 
occurs, it is possible that one waterdrop impact may work-harden a metal 
which has properties comparable to those of 1100-0 aluminum to a suf- 
ficient degree that a crack can form in it when a second impact occurs 
at the same point. If this assumption is made, then, for tetrahedral 
fragments, the minimum value that the first number of the counting 

rule, can have is four; this is one more than the minimum number of 
three cracks required to circumscribe a tetrahedral fragment. The 
numbers of the counting rule, EL, must include both the impacts that 
are required for embrittlement of the metal of the test specimen and 
the crack-forming impacts required to circumscribe a piece of the solid 
and release it as an eroded fragment. 

E. Equation for Rate of Volume Loss 

Rate of volume loss, R, expressed in volume loss per impact is 
R = (volume of material lost) /(number of impacts sustained). 
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If the volumes of eroded fragments are all equal, that is, if the average 
volume of an eroded fragment is used. 


£ _ (number of fragments ejected) (average volume per fragment) 

(number of impacts sustained) 

The quotient of number of fragments ejected divided by the number of 
impacts sustained is the rate of fragment ejection, J. Designating the 
average volume of the typical ejected fragment with the symbol v, the 
rate of volume loss, R, is given by the product 

R = Jv . (4) 

To evaluate R, it is necessary to obtain expressions for J and v. 

1. Rate of Fragment Ejection, J 

Ejection rate, J, is determined both by a statistical function, 

P(x,n) and by a counting rule, . The statistical function gives 
the probability that some number of impacts, x, will occur on the cell 
area, a, when n impacts have occurred on the impingement area A. The 
counting rule, , dictates the number of impacts required to eject a 
fragment from the cell area, a, after some number of fragments (including 
zero) have already been ejected from this particular area. 

a. Function P(x,n) 

The function P(x,n) is uniquely given by binomial probabilities 
if the rain of drops is random. If events are strictly independent in 
the sense that the probability of occurrence of one event does not affect 
the probability of occurrence of a succeeding event, their occurrence is 
described by the well known binomial distribution function. The impacts 
of the drops of a random rain are independent in this sense; the fact 
that a given site has been hit by one drop does not affect the probability 
that it will or will not be hit by a succeeding drop. Consequently, the 
impacts of a truly random rain against the individual cells in an array 
of cells can be exactly described by binomial probabilities. 

The binomial probability, P(x,n), for x impacts on a cell area, 
a, when n impacts have occurred on impingement area A is given by 

P(x,n) = [(n) ! / (x) ! (n-x)!] p X (l-p) n X (5) 
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where p is the probability of a hit on cell area, a. If the random rain 
of drops is uniform, the probability, p, of a hit against any one "of the 
cells having area a is given by 

P = o/A . (6) 

b . Counting Rule 

The magnitude of the counting-rule numbers is determined by a 
number of variables. For a brittle material, the counting rule is primar- 
ily a function of relative impact velocity, V, drop mass, M, angle of 
attack, 0, and the energy per unit volume that the solid is able to 

accept without fracture through the condition that the energy delivered 

1 2 

by a single drop impact, which is given by M (V sin 0) , must exceed 
the energy per unit volume that the solid can accept before fracture will 
result. Clearly, for impacts against a given brittle solid, the magni- 
tude of the counting- rule numbers will be reduced as impact velocity, 
drop mass, and angle of attack are increased. For a given impact velocity, 
drop mass, and angle of attack, the magnitude of the counting-rule numbers 
will be increased as brittle materials of progressively higher strength 
are considered. 

A ductile metal must be work-hardened to the point of embrittle- 
ment before cracks will form in it. The counting-rule numbers for a 
ductile metal include both the impacts required to embrittle it and the 
impacts that are required to circumscribe a fragment with cracks. The 
magnitude of the counting-rule numbers will be larger as the number of 
impacts required to work-harden the metal under test is increased but 
this required number will depend upon the impact velocity, drop mass, 
and angle of attack. 

The magnitude of the counting- rule numbers is also affected by 
conditions that develop as erosion progresses. The counting-rule numbers 
increase in magnitude if dissipation of impact energy occurs as a result 
of drop liquid that remains trapped in the surface roughness. The magni- 
tude of the counting- rule numbers is decreased by the presence of re- 
sidual crack ends and subsurface damage in underlayer material because 
the presence of crack ends and subsurface damage constitute a partial 
loosening of sections of the solid which will eventually break away as 
eroded fragments. 
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c. Function F(x,n) 

The number of fragments ejected from the area A, F(x,n), is 
found by summing the probability, P(x,n) between the limits imposed by 
the counting rule, that is, 

Np + Nj - 1 

F(x,n) = nn T ^ ! P(x,n) + 2rm I 

x = N 0 

N + N, + N„ + No - 1 
o 1 2. £ 

3rm' ^ ) P(x,n) + etc... (7) 

x = N q + N x + N 2 

where P(x,n) is given by Equation 5, and the probability of a hit, p, is 
given by Equation 6. For the necessity of the factor nn', see Section 
•II. B. 3. 

The function F(x,n), when plotted against x or n, is a step 
function because of the limits on the sums which are imposed by the 
counting rule. Let us consider the first term of the function F(x,n) . 

The summation does not start with x = 0, but with x = N Q . No fragment 
can be ejected until at least N c impacts have been received on the area, 
a , of the typical cell; these impacts are received at normal incidence 
(or at some initially prescribed angle of incidence) . 

After these N Q impacts have been received at normal incidence, 
the ejection of the first layer of fragments becomes an accomplished fact. 
This is signified in Figure 5 by the abrupt rise from zero fragments to 
qr)' fragments after N 0 impacts have been received on the cell area, a. 

This abrupt rise is the first riser of the staircase plot of Figure 5. 

At this point a second period of quiescence occurs at impacts from N Q 
to N 0 + N X - 1 are received on the cell area, a. This period of qui- 
escence constitutes the first tread in the staircase plot of Figure 5. 

After the loss of nn* fragments is an accomplished fact, the 
exposed underlayer surface of the cell area, a, is no longer planar; the 
impingement angle is now reduced below 90 degrees. If it is reduced to 
an average value of 45 degrees, the energy delivered per impact is re- 
duced by a factor of two and, on the basis of the reasoning presented 
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N 0 + Ni + N 2 - 1 
Z • P(x,n) + 

x = N q + 
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Number of Fragments Ejected from Area 


4nn ' 


3nn ' 



Number of Impacts in units of N q 


Figure 5. Schematic Representation of the Effect of the Counting Rule 
on the Ejection of Fragments. 
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in Section II. C, N-^ i 2N Q . During the period of quiescence represented 
by the first tread in the staircase plot of Figure 5, impacts are being 
received at an average angle of attack of 45 degrees. 

We are now ready to consider the second term of the function 
F(x,n). After N Q impacts have been received on area a at an angle of 
90 degrees and impacts have been received at the average impingement 

angle of 45 degrees, the ejection of the second layer of fragments from 
the area A becomes an accomplished fact. This is indicated by the second 
riser of the staircase plot of Figure 5. At this point a third period 
of quiescence sets in as the impacts from N 0 + to N 0 + + N 2 - 1 

are received on the cell area, a. This period of quiescence constitutes 
the second tread of the staircase plot of Figure 5. 

After the loss of the second layer of fragments is an ac- 
complished fact, empirical evidence (See Section III) suggests that the 
angle of attack is reduced to a value of about 30 degrees and that the 
counting- rule number N 2 is 4N Q . During the period of quiescence repre- 
sented by the second tread in the staircase plot of Figure 5, impacts 
are being received at an angle of attack of about 30 degrees or with a 
reduced energy that is equivalent to an effective angle of attack of 
30 degrees. 

We are now ready to consider the third term of the function 
F(x,n). After N 0 impacts have been received on area a at an angle of 
90 degrees, impacts have been received at an average impingement angle 
of 45 degrees, and N 2 impacts have been received at an effective impinge- 
ment angle of 30 degrees, the ejection of the third layer of fragments 
from the area A becomes an accomplished fact. This is indicated by the 
third riser in the staircase plot of Figure 5. At this point a fourth 
period of quiescence sets in as impacts from N Q + + N 2 to N 0 +Ni+ 

n 2 + n 3 - 1 are received on the cell area, a. This period of quiescence 
constitutes the third tread of the staircase plot of Figure 5. 

Whether or not the steps of the staircase plot shown in Figure 
5 will be sharply delineated is determined by the magnitude of the proba- 
bility of a hit, p, in conjunction with the magnitude of the counting- 
rule numbers, N^. For a probability, p, of a given size, there is a 
threshold size range for the counting-rule numbers for which the proba- 
bility distribution can be contained on a single tread of the counting-rule 
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staircase. If the size range of the counting-rule numbers is equal to 
or greater than this threshold size, the counting- rule steps will be 
sharply delineated in the plot of fragments lost against number of im- 
pacts sustained. 

If the size range of the counting-rule numbers is smaller than 
the threshold size, the probability distribution will extend over more 
than one tread of the counting-rule staircase. For this condition, the 
counting-rule steps will appear only as oscillations in a plot of 
fragments lost versus number of impacts sustained. If the size range 
of the counting-rule numbers is so small that the probability distri- 
bution extends over many treads of the counting-rule staircase, the step 
nature of the counting rule will not be apparent at all in a plot of 
fragments lost against number of impacts sustained. 

The magnitude of the number of impacts sustained, n, also has 
an effect on the sharpness with which the counting-rule staircase is re- 
flected in the plot of fragments lost against number of impacts sustained. 
This follows because the probability distribution spreads as the number 
of impacts, n, is increased. For values of n that are large enough to 
permit the spread-out probability distribution to occupy more than one 
tread of the counting-rule staircase, the steps of the counting-rule 
staircase will appear only as oscillations in a plot of fragments lost 
against number of impacts sustained; for still larger values of n, the 
steps of the counting-rule staircase will not appear at all in such a 
plot. On the basis of these considerations, a plot of fragments lost 
against number of impacts sustained may exhibit a step structure for 
small values of n and be without this structure for large values of n. 

d. Ejection Rate, Volume Rate, Accumulated Volume Loss 

The average rate of fragment ejection over an interval from 
n^ to n 2 impacts against area A is given by 

J = (F 2 - F x ) / (n 2 - n x ) . (8) 

To make the ejection rate, J, as instantaneous as possible, the quantity 
(n£ - n^) should be made as small as feasible. From Equation 4, the 
average rate of volume loss, R, for each interval of n is Jv where v is 
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the volume of the average ejected fragment (See Section II .E. 2) and J is 
the average rate of fragment ejection. 

The accumulated volume loss can be calculated from the rate of 
volume loss. The increment of volume loss, Q, for each interval of n 
is the product of the rate of volume loss multiplied by the difference 
in n, that is, 

Q = R An . (9) 

The accumulated volume loss, L, up to any value of n is the sum of the 
values of Q up to the specified value of n, that is. 


n 



i = 1 


To calculate the rate of volume loss and the accumulated volume 
loss, it is necessary to have an expression for v, the volume of an 
ejected fragment. 

2. Volume of an Ejected Fragment 

Radial tensile stresses produced by the localized impact pressure 
and radial flow of drop liquid are the means by which an impinging drop 
damages a solid. (1^,16) £ s reasonable to expect that the rate at 

which drop-impact erosion progresses is a function of the variables which 
determine the magnitudes of these quantities. The variables in question 
are the relative impact velocity and the size and mass of the drops that 
impinge. The effect of these variables on the ejection rate, J, was 
discussed in the preceding section. Because relative impact velocity 
and drop mass control the energy delivered per impact, and because the 
energy delivered per impact governs the size of the surface area of the 
cracks that form, the relative impact velocity and drop mass must also 
determine the volume, v, of an eroded fragment. 

a. Nature of the Typical Cell 

( 8 ) 

The initial form of the cellular model v ' was directed to a 
homogeneous material that had uniform bond strength. Such a material is 
ideal; it has no preferred directions for crack propagation and no points 
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of special weakness. The closest approach to an ideal material of this 
kind is a single crystal of high quality. However, even single crystals 
contain defects. The majority of crystals are built up of micromosaic 
units that are 10 - ^ to 10~5 cm in linear dimensions . These units 

are mutually disoriented by angles of the order of seconds or minutes of 

aT .„ ( 17 ) 
arc. 

The boundaries separating the slightly disoriented crystalline 
blocks have been considered to be the more stable ensembles into which 
dislocations have grouped themselves as the result of thermal treatment 
of a deformed crystal. ^9) Bragg and Burgers ^1) W ere the first to 
suggest that the mosaic blocks are units whose bounding surfaces are 
dislocation walls. X-ray and optical investigations have revealed that, 
in addition to the very small micromosaic units, single crystals also 
contain larger macromosaic units which differ slightly in orientation 
and have dimensions of the order of 0.1 to 1 mm and larger. 

The material in dislocation walls around mosaic units may be 
associated with the crack initiation and crack propagation that are 
thought to result in drop- impact erosion. In addition, the material in 
such walls is clearly in a higher energy state than material within the 
mosaic units themselves and, therefore, is more soluble. Chemical 
attack has never been eliminated as a drop-impact-erosion mechanism; 
under impact conditions, even water is a corrosive liquid (See Appendix A). 
In the light of these considerations, it seems reasonable to suppose that 
the typical cell, which becomes an eroded fragment, could be identified 
with the mosaic units of crystalline materials either on the basis of a 
crack-formation mechanism of erosion, a chemical attack mechanism of 
erosion, or a combination of these mechanisms. 

b. Sizes of Typical Eroded Fragments 

The question of the size of cavitation-eroded fragments has 

( 22 ) ( 23 ) 

been explored by Brandenberger and DeHaller , Hoff and Langbein , 

( 9 ) 

and Gould . In general, it was found that most of the eroded fragments 

-4 

have a size very close to that of a micromosaic unit, that is, 10 cm, 
but fragments up to 500 times this size were reported. 

If an eroded fragment is identified with a cr-cell, this evidence 
suggests that it may be possible to associate the micromosaic unit with 
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the typical cell of smallest possible size. That is, the micromosaic 
unit size may constitute the lower limit of the possible size range for 
a-cells . 

It is noteworthy that cavitation-eroded fragments would be 

expected to be smaller than drop- impact-eroded fragments because the 

diameter of the jet formed during the collapse of a cavitation bubble 

is small in comparison with the drop diameters encountered in most drop- 

impact-erosion applications. Also, in this connection, it was observed 
(24) 

by DeHaller that cavitation erosion has a finer texture than drop- 
impact erosion. 

Tentatively, it might be expected that a natural upper limit 
of ejected fragment size may be the grain size of a metal that is under 
test because A.S.T.M. grain sizes from 7 to 0 are, respectively, 30 to 
360 times the 1-micron micromosaic unit size. However, under certain 
circumstances, some eroded fragments may consist of clusters of grains 
and grain fragments. 

If eroded fragments produced in room-temperature tests are 
found to be in the size range of the grains of the test metal, this may 
provide evidence that chemical attack is playing a role in the mechanism 
of failure. Grain boundaries vary in composition from the material with- 
in the grains and are attacked in a corrosive environment. Water, under 
impact conditions, can become a corrosive liquid (See Appendix A). 

(9) 

Gould has stated that the number of large eroded particles 
increases with increasing depth of erosion or with elapsed test time. 

If cell size is identified with eroded-particle size, this suggests 
that, for a calculation of erosion rate that extends over a considerable 
period of time, the typical cell increases in size. The consequences of 
an increase in cell size in terms of a statistical model of erosion rate 
are pointed out briefly in Section IV. 

c. Shapes of Typical Eroded Fragments 

The pyramid having the smallest number of sides in a tetrahedron. 

(8) 

In terms of the initial cellular model of erosion rate'' ' , the eroded 
fragment produced by the least number of hits is a tetrahedral piece of 
the solid which has been cut loose by the intersection of three cracks. 
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Baker, Jolliffe, and Pearson have reported evidence of subsurface 
cracks formed during drop-impact test. If subsurface cracks exist, then 
prismatic fragments or platelets can he cut loose by the intersection of 
three or more cracks that originate at the surface and extend down to 
the subsurface crack lying parallel to the surface. 

2 

An electron micrograph of cavitation-eroded fragments from 
12 percent chromium stainless steel is shown in Figure 6. This micro- 
graph reveals fragments of many different sizes and shapes with a mini- 
mum face diameter of about one micron. The most predominant two- 
dimensional shape appears to be a polygon of at least five sides. As 
seen in two dimensions , some of the fragments are polygons of so many 
sides that they appear to be circular. 

d. Volume of a Tetrahedral Fragment 

The pyramid having the smallest number of sides is a tetra- 
hedron. Let the typical fragment be a tetrahedron circumscribed by 
three semielliptic cracks AB, BC, and CA (See Figure 7) that have origi- 
nated at the surface and have run down into the solid and intersected at 
point D. For simplicity, let the typical fragment be a regular tetra- 
hedron with equal side length. The area of the triangular base, ABC, 

(25) 

of side length b is b^ ^3 74. The altitude h of a regular tetrahedron 
is given by b V"273 and the volume is 1/3 the product of the base 
multiplied by the altitude. Associating the edge length b with the 
minimal or critical tetrahedral fragment having volume, v , 

v c = b 3 y/~P/12 . (11) 

If the edge length b of the minimal tetrahedron is associated with the 

-12 3 

1-micron edge length of a micromosaic unit, then v c = 0.118 x 10 cm . 

The area of one of the intersecting semielliptic cracks that 
form the tetrahedron ABCD in Figure 7 is Trh'b/2 where h' is also the 
slant height of the tetrahedron. Because the slant height of a regular 
tetrahedron is given^-^by b yJV/2, the area of the semielliptic crack 


’This micrograph was obtained by Dr. G. Gould, General Electric Company, 
Materials and Processes Laboratory, Schenectady, New York. 
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Figure 6. Electron Micrograph of Eroded Fragments of 12 Percent 
Chromium Stainless Steel (Photo by Dr. G. Gould). 
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Figure 7. Schematic Drawing of a Tetrahedral Fragment Formed by the 
Intersection of Three Semi-Elliptic Cracks. 
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is irb v? /4. For crack formation in a specific material, if the part 
of the impact energy that is invested in crack formation is doubled, 
the area of each of the intersecting semielliptic cracks will be doubled 
because the surface energy of a crack is the product of the area of the 
crack multiplied by the surface tension of the solid. To double the 
area of each of the semielliptic cracks, the original edge length b must 
increase to V? (b). 

Let E be the part of the impact energy that is invested in 
crack formation and let E c be the critical threshold energy needed to 
form a crack. By analogy with the case in which the area of the semi- 
elliptic cracks is doubled, if the part of the impact energy that is 
invested in crack formation is increased E/E c times, the area of the 
semielliptic cracks is given by (^E/E c Cb) 2 / 4 and the original 

edge length b must increase to yj E/E c b. 

The volume of the tetrahedron cut out by the three larger- 
sized cracks will then be (-^ E/E c (b) 3 y[P/12. If the edge length b is 
associated with the minimal unit of crack length for cracks extending 
from the surface, defined by the critical energy E c or critical velocity 
V c at which a crack is just able to form in the work-hardened state of 
a metal under drop-impact test, then the volume of a typical fragment 
formed at an impact energy associated with the crack formation energy E, 
or impact velocity which corresponds to the crack formation energy E, is 
proportional to (E/E c ) or to (V/V c ) times the volume of the minimal 
tetrahedral fragment associated with the energy E c or with the velocity 

v c . 

For a specific material, the volume of the tetrahedral fragment, 
v t , formed at an impact energy associated with the crack formation energy 
E or at a corresponding impact velocity is 

v t = kv c (V/V c ) 3 (12) 

where k is a proportionality constant and v c is given by Equation 11. 
Associating the edge length b of the minimal tetrahedral fragment with 
the edge length of a 1-micron micromosaic unit, and using Equation 11 
for the volume of the minimal tetrahedral fragment, 

v t = k (0.118 x 10~ 12 ) (V/V c ) 3 . (13) 
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It seems reasonable to expect that the constant k is a function 
of the drop mass, M, as well as of material properties such as the 
strength, S, and the elastic modulus, Y, of the solid being tested; this 
constant must, therefore, be different for different materials, drop 
sizes, and drop-liquid densities. A quantity that might reasonably be 
expected to enter into the expression for the constant k is the tensile 
force needed to open a crack. 


The Griffith force 


(26) 


tending to pull a crack open is given 


by 


and 


1/2 

(2TY/irc) for plane stress 

(2TY/[ttc( 1 - v 2 )] for plane strain 


where T is surface tension, Y is Young's modulus, c is half the length 
of the crack, and v is Poisson's ratio. Because the dimensions of force 
are g.cm/sec^ and because the dimensions of surface tension are g/sec^, 
a dimensionless factor would be obtained if the force tending to pull a 
crack open were to be divided into the product of surface tension multi- 
plied by the drop diameter. 


It seems unlikely that the thermodynamic surface tension should 
be used for T in the Griffith equation because the formation and propa- 
gation of a crack involves more than the creation of new surface. What 
should be used is the surface energy measured in a fracture process; for 
ductile materials such a measurement will be difficult to make. This 
surface energy can be obtained from the strain energy release rate 
(Irwin G-factor) which is proportional to the square of the yield 

strength divided by Young's modulus. For the case of plane stress, if 

2 

T is taken to be proportional to Sy /Y, where Sy is the 0.2 percent off- 
set yield strength of the solid, the Griffith force tending to pull a 
crack open would reduce to a constant times the quotient of yield 
strength divided by the square root of the half-length of the crack. 

If Equation 12 or Equation 13 is substituted into Equation 4, 
it is seen that the volume rate of erosion for a specific material depends 
on the third power of the impact velocity for the case that the typical 
eroded fragment is a tetrahedron. The slant height of a regular pyramid 
of any number of sides is given by [_a^ - (b/2)^J where _a and are 
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the edge lengths of a triangular face. For the case that a = gb wjiiere 

2 1/2 

g is a numerical coefficient, the slant height is b(4g - 1) /2. Be- 

cause the slant height of a regular pyramid of any number of sides is 
given by the same power of the edge length b as the slant height of a 
regular tetrahedron, it appears that the same velocity dependence of 
the rate of erosion should be obtained for pyramidal fragments of any 
number of sides as is obtained for tetrahedral fragments. 

e. Volume of a Prismatic Fragment 

The existence of subsurface cracks in eroded specimens was 
reported by Baker, Jolliffe, and Pearson^"^ ; however, the number of 
such cracks appeared to be restricted. If subsurface cracks that run 
parallel to the impacted surface are produced either as a result of 
subsurface shear stresses or by the growth of cracks from the side walls 
of existing craters, then cracks that run down from the surface may 
intersect the subsurface cracks to form eroded fragments that are prisms 

The potential energy invested in a crack of fixed depth is 
proportional to the length of the crack. Assuming that the fraction, E, 
of the impact energy that is used in crack formation remains constant 
as the impact energy is increased and that the depth of penetration of 
a crack extending downward from the surface is fixed by a subsurface 
crack running parallel to the surface, then, if the impact energy is 
doubled, the length of the crack which extends downward from the surface 
to the subsurface crack should also be doubled. 

Let the typical fragment be a triangular prism circumscribed 
by three cracks, a, b, and c, and let it have an altitude or thickness 
determined by a pre-existing subsurface crack (See Figure 8, upper view) 
The area of the triangular upper and lower base, ABC, is given by 
£s (s - a) (s - b) (s - c)] where s = (a + b + c)/2. If the impact 
energy is doubled, then the side lengths a, b, and c, will be doubled 
and the area of the upper and lower base will be increased by a factor 

p 

of 2* or 4 times. 

Similarly, for the case that the fragment is a polyhedron of 
four sides (See Figure 8, lower view) also having altitude or thickness 
determined by a pre-existing subsurface crack, if the impact energy is 
doubled then the lengths of the sides a, b, c, and d are doubled and 
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the area circumscribed by the four cracks is increased by a factor of 2 
or 4 times. In general, if the impact energy is increased z-fold, the 
area circumscribed by the cracks that result is increased z^ times. 

If a minimal unit of crack length extending from the surface 
to the subsurface crack is defined by the critical energy, E c , or critical 
velocity, V c , at which such a crack is just able to form, then the surface 
area of a typical fragment formed at an energy, E, or at a velocity V, is 

o / 

proportional to (E/E c ) or to (V/V times the minimal area associated 
with the energy, E c , or with the velocity, V c . The minimal area is the 
area of a fragment of the smallest possible size. This minimal area is 
designated as A c . 

In view of the preceding considerations, it appears that the 
volume of a prismatic eroded fragment, Vp, should be proportional to 
(A c ) (t 0 ) (V/V c ) ^ where t CT is the depth of the subsurface crack below the 
surface, that is 

v p = k'A c t CT (V/V c ) 4 , (14) 

where k' is a proportionality constant which is a function of drop mass 
and material properties. If Equation 14 is substituted into Equation 4, 
it appears that the volume rate of erosion for a specific material de- 
pends on the fourth power of the impact velocity for the case that the 
eroded fragments are prismatic. Because the area of all regular polygons 
is a function of the square of the edge length and because the areas of 
irregular polygons can be expressed as the sums of the areas of regular 
polygons, it appears that the same velocity dependence of the rate of 
erosion should be obtained for prismatic fragments of all sizes and 
shapes . 


III. PARTIAL TEST OF THE EQUATIONS 

Equations 4, 7, 13, and 14 are not ready for independent test be- 
cause they contain quantities that cannot yet be evaluated from theoreti- 
cal considerations alone. They can be given a partial test by assessing 
the values of the needed quantities with recourse to experimental data. 

It was anticipated that partial test of these equations would in itself 
provide valuable information. 
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Both drop-impact and cavitation weight-loss data were collected 

under very carefully controlled conditions for the purpose of testing 

equations. The test specimens were of zinc, iron, tantalum, nickel, 

Udimet 700 in the solutioned state and Udimet 700 in the aged state. 

These metals represent the face-centered-cubic, body-centered-cubic, 

and hexagonal lattice-packing types. Information on the work-hardening 

characteristics, tensile stress-strain characteristics, impact behavior, 

and minimum velocity required to dent by impact of a steel sphere was 

(14) 

obtained for these metals in an earlier study. The specimens of 

these metals that were used for drop- impact and cavitation erosion tests 

were from the same heat and were in the same heat-treatment state as the 

(14) 

impact specimens used in the earlier study. 

The experimental values of weight loss for these metals were de- 
termined after arbitrary periods of waterdrop impingement on a rotating- 
3 

arm device. All of the tests were performed at a relative impact ve- 
locity of 1000 ft/sec with 0.0866-cm-diameter drops. With the apparatus 
that was used, 1 both the impingement velocity and the drop diameter 
could be held constant with a high degree of precision. With the ex- 
ception of Udimet 700, which was found to require long periods of test 
before erosion loss occurred, a separate test specimen was used for each 
data point. This eliminated changes in the impingement area. A, which 
occur as a consequence of removal and replacement of a test specimen in 
order to determine weight loss. 

The drops impinged within a more or less circular area having a 
radius equal to the drop diameter. The reported rate at which the drops 
impinged was 161.5 per second. The reported data for the tests performed 
were test time, initial weight of specimen, and final weight of specimen. 
For a detailed description of the test apparatus, a discussion of pro- 
cedures used in carrying out the tests, and a tabulation of the raw data, 
see Appendix B. 

Weight loss during each time of exposure to drop impingement was 
found by subtracting the weight of the specimen after exposure from its 


The tests, which were obtained from the University of Minnesota, were 
carried out under the direction of Mr. John A. Almo; Mr. Michael McKay 
carried out most of the test runs. 
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weight before exposure. Weight loss was divided by the density of the 
metal to convert it to volume loss. Test duration times were multiplied 
by the rate at which the drops impinged to convert them to numbers of 
drop impacts sustained. The resulting data for the six metals are given 
in Table 1. 

It was not possible to include a test of the equations with use of 
the cavitation data, which were collected independently by Mr. Stanley 
G. Young at the NASA Lewis Research Center, Cleveland, Ohio, in the 
present study. This test and whatever correlation may be found between 
the drop-impact and cavitation erosion data that were collected will be 
presented in a later study. 

A. Relationships Needed to Test Equations 4 and 7 

In the development of the statistical treatment, it was assumed 
that the typical cell in the array of cells has the size of an eroded 
fragment and that the diameter, d, of a 0 -cell is equal to the distance 
across an eroded fragment. From Equation 6, taking the area of a o-cell 
in the direction of impingement to be circular, 

p = d 2 7r/4A 

from which 

d = (4A p/ it) 1/2 . (15) 

Let it now be assumed that the number of impacts required to work- 
harden a piece of the solid, circumscribe it with cracks, and release 
it as an eroded fragment is proportional to the size of the fragment. 

The number of impacts required to release a fragment from the original 
surface of the test specimen at the start of a test is N 0 . Then N Q is 
proportional to d. By trial, it appears from recourse to the experi- 
mental data that the proportionality constant is the reciprocal of the 
edge length of a micromosaic unit or about 10^ cm - -*-. This leads to the 
relationship , 

N = 10 4 d (16) 

o 

which states that the number of impacts required to release a fragment 
is equal to the diameter of the fragment expressed in microns. 
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TABLE 1 


EXPERIMENTAL VOLUME LOSS DATA FOR THE SIX SELECTED METALS 


ZINC 

NICKEL 

IRON 

TANTALUM 

UDIMET 700 

Sol'n. 

UDIMET 700 Aged 

Number of 

Volume 

Number of 

Volume 

Number of 

Volume 

Number of 

Volume 

Number of 

Volume 

Number of 

Volume 

Impacts 

Loss 

Impacts 

Loss 

Impacts 

Loss 

Impacts 

Loss 

Impacts 

Loss 

Impacts 

Loss 

, „3 t 

10 Impacts 

-4 

10 cc 

4 

10 Impacts 

-4 

10 cc 

5 

10 Impacts 

-4 

10 cc 

10^ Impacts 

-4 

10 cc 

0 

10 Impacts 

-4 

10 cc 

6 

10 Impacts 

-4 

10 cc 

2.42 

0.70 

0.236 

1.01 

0.388 

1.77 

0.563 

0.00 

1.464 

3.03 

1.454 

0.37 

2.46 

0.42 

0.491 

0.22 

0.775 

3.04 

0.969 

2.77 

1.755 

7.82 

2.038 

6.06 

4.84 

1.40 

0.581 

1.46 

0.969 

6.10 

1.292 

3.25 

2.046 

9.21 

2.328 

9.72 

4.87 

1.54 

1.000 

2.69 

1.132 

6.99 

1.948 

5.96 

2.627 

12.50 

2.910 

13.38 

9.69 

2.80 

1.455 

1.46 

1.295 

7.24 

2.422 

4.64 

3.208 

14.39 

3.491 

16.41 

9.74 

2.24 

2.426 

3.48 

1.551 

6.98 

3.230 

8.37 

4.371 

18.30 

4.072 

18.93 

12.23 

5.75 

4.882 

7.07 

1.950 

7.24 

3.908 

12.53 

5.419 

23.61 

4.654 

21.46 

14.54 

5.60 

9.693 

12.8 

2.326 

9.52 

5.249 

15.54 

6.687 

29.16 

5.235 

23.23 

14.54 

5.74 

14.535 

16.62 

3.101 

11.68 

7.751 

17.29 

8.011 

33.45 

5.814 

25.00 

19.38 

6.30 

19.383 

15.61 

4.652 

15.74 

9.690 

19.04 

9.303 

34.97? 

6.977 

28.91 

19.44 

6.02 

24.225 

22.12 

6.202 

18.41 

11.144 

29.70? 


hole 

8.140 

30.68 

21.93 

15.42 

29.070 

19.43 

8.075 

26.42 

11.306 

21.45 



9.302 

33.08 

24.23 

15.70 

33.918 

23.14 

9.690 

27.94 

15.504 

25.42? 



10.465 

35.35? 

24.63 

15.42 

38.76 

23.14 




hole 




hole 

29.07 

16.54 

43.62 

30.55? 









29.07 

17.38 











29.51 

25.51? 











31.51 

22.01 











33.92 

21.02 











34.00 

26.49 











36.34 

30.56 











36.35 

29.72 











38.76 

35.04 











38.79 

34.90 











41.18 

34.34 











41.52 

32.52? 












hole 
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Substituting Equation 15 into Equation 16, 

N q = 1.12838 x 10 4 (Ap) 1/2 . (17) 

For the experimental data that were collected, A = 0.0235606 cm^. For 
these data only . 

1/2 

N q = 1732 p ' (Rl) 

and 

p = (N 0 /1732) 2 . (R2) 

The numbers given to these equations indicate that they are restricted 
in application to an impingement area of 0.0235606 cm^. 

The number of cells of area a in an array that extends over area A 
is given by Equation 1. If a rain of n impacts occurs over area A, then 
the average number of impacts against a typical a-cell in the array is 
the quotient na/A. With use of Equation 6, the average number of impacts 
that occur against the typical cell when a rain of n impacts occurs over 
area A is just p n where p is the probability of a hit against the cell. 

A layer of nn' fragments will be released from the area A when the 
number of impacts that occur at a typical cell site is N Q or when the 
average number of impacts is p n Q and n Q is a specific number of impacts 
that occur over area A. Remembering this fact, and multiplying Equation 
R2 by n Q , 

pn = N 2 n / (1732) 2 = N . 

* O O O O 

From the second equality, 

n Q = (1732) 2 /N q = 2999824/N q . (R3) 

B. Test of the Equations with Experimental Drop-Impact Volume-Loss 
Data for Iron 

The drop-impact volume-loss data that were collected for iron are 
particularly well suited for test of the equations because a large number 
of data points happened to be collected over the range of impacts in which 
the maximum rate of volume loss occurred. These experimental volume-loss 
data for iron are plotted against number of impacts sustained in Figure 9. 
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1. Expected Form of the Rate-of-Volume-Loss Curve 

The maximum rate of volume loss occurs at the number of impacts for 
which the first sharp rise in accumulated volume loss takes place. In 
Figure 9, it can be seen by inspection that this first sharp rise in 
accumulated volume loss takes place somewhere between 0.9 and 1 x 10^ 
impacts. Further inspection of Figure 9 indicates that two other sharp 
rises in accumulated volume loss also occur. These rises, which are 
located at about 2.75 x 10^ and 6.5 x 10^ impacts, respectively, are 
associated with peaks in the rate of volume loss that are lower than the 
first peak; the first peak is the maximum rate of volume loss. 

The peaks in the rate of volume loss plotted against impacts sus- 
tained are separated by regions in which essentially no volume loss 
occurs. Three of such flat plateaus have been roughly indicated in 
Figure 9. The rate of volume loss on such a plateau is zero. From 
inspection alone, it can be concluded that the rate of volume loss for 
iron plotted against number of impacts sustained will have the form 
indicated schematically in Figure 10. 

For a further discussion of features of rate curves see Section III.G. 

2. Sample Calculation of Rate-of-Volume-Loss and Accumulated- 
Volume-Loss Curves 

The theoretical work has not yet been carried far enough to pro- 
vide an expression for the first number of the counting rule, N Q , or 
for the number of impacts, n Q , that must occur against an area. A, of 
a given material before the first sharp rise in accumulated volume loss 
takes place. It is necessary to have recourse to the experimental data 
to obtain the number of impacts, n Q . On the basis of the experimental 
data, 0.95 x 10^ can be accepted as a probable value of the impacts 
that are required to reach the first sharp rise in accumulated volume 
loss for iron. Substituting this value of n Q into Equation R3, the 
value of N q is found to be 32. Substituting this value of N 0 into 
Equation R2, the probability of a hit, p, is found to be 3.4135 x 10 - ^. 

For iron, with N 0 equal to 32, the counting rule numbers, N^, can be 
expected to pass through the values = 64 and ^ = 128 as the angle 
at which the drops impinge changes from normal incidence to an effective 
angle of 30 degrees (See Section II. C). 
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Rate of Volume Loss, cc/ impact 
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The input data for iron to be used in the computer program (See 

--4 

Appendix C-B) for Equation 7 are: p = 3.4135 x 10 , N Q = 32, 

Nl = 64, and N 2 = 128. Because it was desired to calculate the rate of 
erosion for iron, a printout of the accumulated number of ejected layers 
of fragments was required and the condition = R 2 = R 3 “ 1 was imposed 
(See Appendix C-B) . The computer printout of layers of fragments 
ejected is given in the second column of Table 2. 

From Equation 16 and the value of 32 for N Q , the diameter, d, of 
an eroded iron fragment is 0.0032 cm or 32 microns. This is about half 
the average grain size of the iron metal used, which was found to be 
63 microns using the Heyn intercept method (ASTM Designation E112-63) . 
Assuming that the area of a typical a-cell in iron is the area of a 
circle having the diameter d = 0.0032 cm, the area, a, of the typical 
cell or eroded fragment for iron is 8.0425 x 10 -D cm . Remembering that 
the impingement area. A, for the drop-impact tests performed is 0.0235606 

O 

cm , the number of cr-cells, q, for iron from Equation 1 is 2929.51. For 
the 0. 0866-cm-diameter drops that were used, n' for iron from Equation 2 
is 42.5098, and the product qq' is 124532. The fragments ejected from 
area A when each a-cell in this area has received some integral number 
of impacts is the product of nn * times the number of layers ejected. 

The fragments ejected from area A are listed in column three of Table 2 
and the values of average ejection rate, J, found with use of Equation 8 
are given in column four. 

The volume of a tetrahedral fragment is 0.11785 times the cube of 
the edge length. If the value d = 0.0032 cm is used for the edge length, 
the volume of a tetrahedral fragment for iron is 3.8617 x 10“® cm^. In 
order to make the first knee of the theoretical volume-loss-versus- 
impacts-sustained curve occur at the same value of accumulated volume 
loss as the first knee formed by the experimental data points, it was 
found by trial that the volume of the typical fragment for iron is 1.5 
times this tetrahedral volume. This suggests that the ejected fragments 
may not all be tetrahedra. It is noteworthy that an eroded fragment 
volume larger than the tetrahedral volume is reasonable because some of 
the eroded fragments may be prisms; an eroded fragment volume smaller 
than the tetrahedral volume could not be justified. With use of the 
factor 1.5, the volume of an average iron fragment is given by the 
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TABLE 2 


SUMMARY OF DATA FOR THE CALCULATED CURVES FOR IRON 


Impacts 

Layers Ejected 

Fragments Ejected 

Ejection Rate 

Volume 

Rate 

Volume Increment 

Accumulated 

Volume 

r 




From Area A 

J 

R 



Q 


Loss , 

L 




— 

— 

frgs/impact 

cc/ impact 

CC 

_ 


CC 


0.25 

X 

10 5 

0.0000 

0 

0.00000 

0.00000 



0.0000 



0.0000 



0.50 

X 

10 5 

0.0008 

99 

0.00025 

0.00143 

X 

10 -9 

0.0001 

X 

10 -4 

0.0001 

X 

10 -4 

0.60 

X 

i0 5 

0.0110 

1374 

0 . 1275 

0.7387 

X 

10 -9 

0.0739 

X 

10 -4 

0.0740 

X 

10 -4 

0.70 

X 

10 5 

0.0648 

8071 

0.6697 

3.879 

X 

10 -9 

0.3879 

X 

10 -4 

0.4619 

X 

1°_4 

0.75 

X 

10 5 

0.1236 

15398 

1.465 

8.488 

X 

10 -9 

0.4244 

X 

10 -4 

0 . 8863 

X 

10 -4 

0.80 

X 

10 5 

0.2078 

25872 

1.780 

10.31 

X 

10 -9 

0.5156 

X 

10 -4 

1.402 

X 

10_4 

0.89 

X 

10 5 

0.4082 

50836 

2.774 

16.07 

X 

10 -9 

1.446 

X 

10 -4 

2.848 

X 

l °-4 

0.90 

X 

10 5 

0.4325 

53864 

3.028 

17.54 

X 

10 -9 

0.1754 

X 

10 -4 

3.023 

X 

10 -4 

0.91 

X 

10 5 

0.4569 

56901 

3.037 

17.59 

X 

10 -9 

0.1759 

X 

10 -4 

3.199 

X 

10 -4 

0.92 

X 

10 5 

0.4813 

59936 

3.035 

17.58 

X 

10 _9 

0.1758 

X 

10 -4 

3.375 

X 

10 -4 

0.93 

X 

10 5 

0.5056 

62958 

3.022 

17.51 

X 

10 -9 

0.1751 

X 

io 4 

3.550 

X 

1°_4 

0.94 

X 

10 5 

0.5 296 

65956 

2.998 

17.37 

X 

10 -9 

0.1737 

X 

io 4 

3.724 

X 

io 4 

0.95 

X 

10 5 

0.5534 

68919 

2.963 

17.16 

X 

10 J 

0.1716 

X 

10 -4 

10 -4 

3.896 

X 

i°4 

0.96 

X 

10 5 

0.5769 

71838 

2.919 

16.91 

X 

10 -9 

0.1691 

X 

4.065 

X 

10 -4 

0.97 

X 

10 5 

0.5999 

74704 

2.866 

16.60 

X 

10 -9 

0.1660 

X 

10 -4 

4.231 

X 

1°_4 

0.98 

X 

10, 

0.6224 

77508 

2.804 

16.24 

X 

10 -9 

0.1624 

X 

10 -4 

4.393 

X 

l0 -4 

0.99 

X 

10 

0.6444 

80242 

2.734 

15.84 

X 

10 -9 

0.1584 

X 

io 4 

4.551 

X 

10 -4 

1.00 

X 

10 5 

0.6657 

82901 

2.659 

15.40 

X 

10 _9 

0.1540 

X 

10 -4 

4.705 

X 

10 4 

1.10 

X 

10 5 

0.8383 

104399 

2.150 

12.45 

X 

10 -9 

1.245 

X 

10 -4 

5.951 

X 

10 _4 

1.20 

X 


0.9351 

116445 

1.205 

6.978 

X 

10 -9 

0.6978 

X 

io 4 

6.649 

X 

10 -4 

1.25 

X 

10 5 

0.9614 

119727 

0.6564 

3.802 

X 

10 -9 

0.1901 

X 

10 -4 

6.839 

X 

10 -4 

1.50 

X 

10 5 

0.9984 

124328 

0.1840 

1.066 

X 

10 -9 

0.2665 

X 

10 -4 

7.105 

X 

1°_ 4 

1.75 

X 

10 5 

1.000 

124529 

0.00804 

0.04657 

X 

10 - 

0.0116 

X 

io 4 

7.117 

X 


2.00 

X 

10 5 

1.001 

124642 

0.00452 

0.02618 

X 

10 -9 

0.0065 

X 

10 -4 

7.123 

X 

10 -4 

2.25 

X 

10 5 

1.019 

126907 

0.09060 

0.5248 

X 

10 -9 

0.1312 

X 

10 -4 

7.254 

X 

1°_ 4 

2.50 

X 

X °5 

1.136 

141503 

0.5838 

3.382 

X 

10 -9 

0.8455 

X 

l °-4 

8.100 

X 

1°_ 4 

2.75 

X 

X °5 

1.427 

177668 

1.447 

8.380 

X 

10 -9 

2.095 

X 

10 -4 

10.19 

X 

1°_ 4 

2.78 

X 

X °5 

1.468 

182871 

1.734 

10.05 

X 

10 -9 

0.3014 

X 

10 -4 

10.50 

X 

10 -4 

2.80 

X 

10 5 

1.496 

186346 

1.738 

10.06 

X 

10 -9 

0.2013 

X 

10 -4 

10.70 

X 

10 -4 

2.85 

X 

10 5 

1.565 

194937 

1.718 

9.953 

X 

10 -9 

0.4976 

X 

10 -4 

11.20 

X 

l °-4 

2.90 

X 

10 5 

1.632 

203192 

1.651 

9.564 

X 

10 -9 

0.4782 

X 

10 -4 

11.67 

X 

1°_ 4 

2.95 

X 

10 5 

1.694 

210899 

1.541 

8.929 

X 

10 -9 

0.4464 

X 

10 -4 

12.12 

X 

10-4 

3.00 

X 

10 

1.750 

217898 

1.400 

8.108 

X 

10 

0.4054 

X 

10 

12.53 

X 

10 










TABLE 2 (Continued) 


SUMMARY OF DATA FOR THE CALCULATED CURVES FOR IRON 


Impacts 

n 

Layers Ejected 

Fragments Ejected 
From Area A 

Ejection Rate 
J 

f rgs/ impact 

Volume Rate 
R 

cc/ impact 

Volume Increment 
Q 
cc 

Accumulated Volume 
Loss , L 
cc 

3.05 x 10g 

1.799 

224087 

1.238 

7.170 x 10_g 

-4 

0.3585 x 10 . 

12.88 x 10"? 

-4 

3.10 x 10° 

1.842 

229419 

1.066 

6.177 x 10 „ 

0.3089 x 10 , 
-4 

13.19 x 10 4 

3.15 x 10 

1.878 

233899 

0.8960 

5.190 x 10_ g 

0.2595 x 10_ 4 

13.45 x 10_ 4 

3.25 x 10 

1.931 

240515 

0.6616 

3.832 

x 10 _ 

0.3832 x 10 , 

13.84 

x 10 . 

3.50 x 10° 

1.988 

247572 

0.2823 

1.635 x 10_ 9 

0.4088 x 10“: 
-4 

14.24 x 10 , 

-4 

3.75 x 10 c 

1.999 

248892 

0.05280 

0.3058 x 10 

0.0765 x 10 . 

-4 

14.32 x 10 4 

4.00 x 10 K 

2.000 

249050 

0.00632 

0.0366 x 10 

0.0092 x 10 . 

-4 

14.33 x 10 

4.25 x 10 c 

2.000 

249063 

0.00052 

0.0030 x 10_ q 

0.0008 x 10_. 

14.33 x 10 ’ 
-4 

4.50 x 10 

2.000 

249063 

0.00000 

0.0000 x 10 9 

0.0000 x io _ 4 

14.33 x 10, 

4.75 x 10g 

2.000 

249064 

0.00004 

0.0002 x 10 . 

0.0001 x io _ 4 

14.33 x 10 , 
-4 

5.00 x 10,. 

2.000 

249070 

0.00024 

0.0014 x 10_ g 

0.0003 x 10 , 
-4 

14.33 x 10 . 

-4 

5.25 x 10° 

2.001 

249150 

0.00320 

0.0185 

x 10 . 

0.0046 x 10 

14.34 

x 10 . 

5.50 x 10° 

2.006 

249744 

0.02376 

0.1376 x 10 g 

0.0344 x 10 ^ 

14.37 x 10 4 

5.75 x 10° 

2.028 

252537 

0.1117 

0.6471 x 10 g 

0.1618 x 10 4 

14.53 x 10 

6.00 x 10° 

2.097 

261149 

0.3445 

1.995 x 10 g 

0.4989 x 10 ^ 

15.03 x 10 _ * 

6.25 x 10° 

2.242 

279142 

0,7197 

4.169 x 10 g 

1.042 x 10~^ 

16.07 x 10 ^ 

6.30 x 10° 

2.280 

283889 

0.9494 

5.499 x 10 g 

0.2750 x 10 * 

16.35 x 10 I 

6.40 x 10° 

2.363 

294256 

1.037 

6.005 x 10 g 

0.6005 x 10 4 

16.95 x 10 4 

6.50 x 10° 

2.452 

305382 

1.113 

6.445 x 10 g 

0.6445 x 10 4 

17.59 x 10 ^ 

6.60 x 10° 

2.543 

316703 

1.132 

6.558 x 10g 

0.6558 x 10 ^ 

18.25 x 10 ^ 

6.70 x 10° 

2.631 

327641 

1.094 

6.336 x 10 g 

0.6336 x 10 4 

18.88 x 10 4 

6.75 x 10° 

2.672 

332803 

1.032 

5.980 x 10 g 

0.2990 x 10 4 

19.18 x 10 ^ 

6.80 x 10° 

2.712 

337690 

0.9774 

5.662 x 10 g 

0.2831 x 10 ^ 

19.46 x 10 

6.90 x 10° 

2.782 

346484 

0.8794 

5.094 x 10 g 

0.5094 x 10 * 

19.97 x 10 ^ 

7.00 x 10° 

2.841 

353822 

0.7338 

4.251 x 10 g 

0.4251 x 10 4 

20.40 x 10 4 

7.25 x 10° 

2.938 

365893 

0.4828 

2.797 x 10 g 

0.6992 x 10 4 

21.10 x 10"^ 

7.50 x 10° 

2.981 

371184 

0.2116 

1.226 x 10 g 

0.3065 x 10 4 

21.40 x 10 * 

7.75 x 10° 

2.995 

372986 

0.07208 

0.4175 x 10 g 

0.1044 x 10 4 

21.51 x 10 ^ 

8.00 x 10° 

2.999 

373470 

0.01936 

0.1121 x 10 g 

0.0280 x 10 * 

21.54 x 10 

8.25 x 10° 

3.000 

373578 

0.00432 

0.0250 x 10 g 

0.0063 x 10 

21.54 x 10 

8.50 x 10? 

3.000 

373622 

0.00176 

0.0102 x 10_g 

0.0025 x 10 - ^ 
-4 

21.55 x 10“; 

-4 

8.75 x 10 

3.001 

373773 

0.00604 

0.0350 x 10 

0.0087 x 10. 

21.55 x 10_ 4 

9.00 x 10° 

3.007 

374416 

0.02572 

0 . 1490 

x 10 „ 

0.0372 x 10 

21.59 

x 10 . 

9.25 x 10° 

3.024 

376529 

0.08452 

0.4896 x 10 g 

0.1224 x 10 4 

21.71 x 10" 4 

9,50 x 10° 

3.067 

381909 

0.2152 

1.247 x 10 g 

0.3116 x 10 4 

22.03 x 10'1 

9.75 x 10° 

3.153 

392638 

0.4292 

2.486 x 10 „ 

0.6215 x 10 , 
-4 

22.65 x 10"? 

-4 

10.00 x 10 

3.289 

409609 

0.6788 

3.932 

x 10 

0.9831 x 10 

23.63 

x 10 





product of 0.18856 times the cube of the edge length. Remembering that 
the volume of a cube is unity times the cube of the edge length, and that 
the volume of a tetrahedron is 0.11785 times the cube of the edge length, 
it appears that the eroded iron fragments are largely tetrahedral. 

The rates of volume loss for iron found with use of Equation 4 and 
with the volume of an eroded fragment, v, equal to 1.5 times the volume 
of tetrahedron with an edge length of 0.0032 cm are given in column five 
of Table 2. The accumulated volume losses can be back-calculated from 
the rates of volume loss to provide a check of the curve of rate-of- 
volume-loss plotted against number of impacts sustained. The increments 
of volume loss, Q, calculated with use of Equation 9 are given in column 
six of Table 2. The accumulated volume losses calculated with use of 
Equation 10 are given in column seven. 

Both the rate of volume loss and the accumulated volume loss found 
from Equation 7 with use of the computer program are shown in Figure 11 
along with the experimental data points. Because the calculated rate 
of volume loss is an average value for an interval of impacts, An, the 
calculated values of rate of volume loss are plotted at the midpoints 
of the intervals. Expressed in time, the interval lengths are from 
6 to 155 sec. It can be seen from Figure 11 that the theoretical volume 
loss back-calculated from the theoretical rate of volume loss is in good 
agreement with the experimental data points up to 7 x 10"* impacts ; for 
larger numbers of impacts sustained, the calculated curve falls below 
the data points. This appears to be the result of the fact that the 
eroded fragment size is increasing (See Section II.E.2.b) whereas the 
theoretical curves have been calculated for eroded fragments of equal 
size. Use of Equation 7 to calculate the accumulated volume loss curve 
for fragments of increasing size is discussed in Section IV. 

C. Test of the Equations with Experimental Drop-Impact Volume-Loss Data 
for Nickel 

By inspection of the experimental data, the first sharp rise in 
accumulated volume loss for nickel was taken to occur at n = 30,000 
impacts. From Equation R3, the value of N Q is 99.9941 and from Equation 
R2 the probability, p, of a hit against the typical cell is 0.00333314. 
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Figure 11. Theoretical Rate-of-Volume-Loss and Accumulated-Volume-Loss Curves for 
Iron Calculated with Fragments of Equal Volume. 


Accumulated Volume Loss, cc 



The counting rule, N^, was expected to pass through the values N-^ = 
198.9882 and N£ = 399.9764 as the angle at which the drops impinge changes 
from normal incidence to 30 degrees. From inspection of the experimental 
data, it was concluded that an angle of attack smaller than 30 degrees 
was reached. Because = 2N Q and N 2 = 4N 0 , the value of N 3 was taken 
to be 8 N Q or 799.9528. 

From Equation 16 and the value of 99.9941 for N 0 , the diameter, d, 
of an eroded nickel fragment was found to be 100 microns. This is roughly 
double the average grain size of the nickel metal used which was found to 
be 53 microns with use of the Heyn intercept method. The volume of a 
tetrahedral fragment is 0.11785 times the cube of the edge length. The 
calculated value of a tetrahedral fragment for nickel is 1.1785 x 10“^ cm^. 
The average volume of an eroded nickel fragment was taken to be 1.49 times 
the tetrahedral volume or 1.755965 x 10 cm-' because this fragment volume 
was required in order to bring the first knee of the theoretical volume 
loss curve into agreement with the first knee formed by the experimental 
data points. 

Both the rate of volume loss and the accumulated volume loss found 
for nickel with use of the computer program for Equation 7 are shown in 
Figure 12 along with the experimental data points. It can be seen from 
Figure 12 that the theoretical volume loss backcalculated from the theo- 
retical rate of volume loss is in relatively good agreement with the ex- 
perimental data points. The cluster of data points near the origin of 
this plot suggests that the first sharp rise in accumulated volume loss 
for nickel should perhaps have been taken to occur at less than 10,000 
impacts. This would make the maximum rate of volume loss for nickel 
larger by nearly a factor of two. However, such a choice would result 
in a value of N Q which is not compatible with the remaining numbers of 
the counting rule. 

D. Test of the Equations with Experimental Drop-Impact Volume-Loss Data 
for Tantalum 

By inspection of the experimental data, the first sharp rise in 
accumulated volume loss for tantalum was taken to occur at about n Q = 
120,000 impacts. On the basis of this value of n Q , the value of the 
first number of the counting rule, N Q , was found to be 24.9986 and the 
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Figure 12 


Theoretical Rate-of-Volume-Loss and Accumulated-Volume-Loss Curves for Nickel 
Calculated with Fragments of Equal Volume. 


Accumulated Volume Loss, cc 



probability, p, of a hit against the typical cell was found to be 
0.000208321. The counting rule, N^, was expected to pass through N-^ = 
49.9970 and N 2 = 99.9940 as the angle at which the drops impinged changed 
from normal incidence to 30 degrees. 

The diameter, d, of an eroded fragment was found to be 24.9986 
microns. This is about 1/5 of the average grain size which was found 
to be 126 microns with use of the Heyn intercept method. The calculated 
volume of a tetrahedral fragment of tantalum was found to be 1.841097 x 

_ Q Q 

10 cm- 3 . It was found necessary to take the average volume of an eroded 

tantalum fragment to be 1.07 times the tetrahedral volume or 1.969974 x 
-9 3 

10 cm in order to bring the first knee of the theoretical volume loss 
curve into agreement with the first knee formed by the experimental data 
points. 

The rate of volume loss and the accumulated volume loss found for 
tantalum with use of the computer program for Equation 7 are shown in 
Figure 13 along with the experimental data points. From Figure 13 it 
can be seen that the agreement between the theoretical volume loss back- 
calculated from the theoretical rate of volume loss is in reasonably 
good agreement with the experimental data points. In the light of re- 
sults reported in Section IV, the relatively poor agreement from about 
4 x 10^ to 8 x 10^ impacts appears to be due to an increase in the volume 
of an eroded fragment. 

E. Test of the Equations with Experimental Drop-Impact Volume-Loss Data 
for Zinc 

The first sharp rise in accumulated volume loss for zinc was taken 
to be n Q = 10,000 impacts. On the basis of this choice, the first number 
of the counting rule, N q , is 300, and the probability, p, of a hit against 
the typical cell is 0.0300018. The diameter of an eroded zinc fragment is 
300 microns. The average grain size of the zinc metal used for the test 
specimen was found to be just 300 microns. From this it appears that the 
zinc eroded in units of entire grains. 

Inspection of the step structure of the experimental accumulated 
volume loss data for zinc indicates that there was no change in the angle 
of attack of the impinging drops, that is, N Q = = N 2 ... = 300. There 

is no obvious answer for the fact that zinc is different from the other 
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Figure 13. Theoretical Rate-of-Volume-Loss and Accumulated-Volume-Loss Curves 
for Tantalum Calculated with Fragments of Equal Volume. 


Accumulated Volume Loss, cc 



metals In this respect. However, the zinc rod was formed by extrusion 
and an X-ray study of the metal showed that it possessed some orientation; 
the hexagons of the lattice structure would be seen if the rod were viewed 
end on. The axis of the test specimen coincided with the axis of the rod 
stock. The fact that whole grains were broken out of the hexagonal 
structure that would be seen if a test specimen were viewed end on may 
account for the fact that there was no change in the angle of attack of 
the impinging drops. 

The rate of volume loss and the accumulated volume loss found for 
zinc with use of the computer program for Equation 5 are shown in Figure 
14 along with the experimental data points. From Figure 14 it can be 
seen that the agreement between the theoretical curve and the experi- 
mental data is good only for the initial section of the plot. In the 
light of results reported in Section IV the poor agreement found for 
numbers of impacts greater than 20,000 appears to be caused by a very 
marked increase in the size of the eroded fragments. 

F. Test of the Equations with Experimental Drop-Impact Volume-Loss Data 
for Aged Udimet 700 Alloy 

Udimet 700 is a high-strength alloy. Failure of such a material 
can be expected to start at weak points (See Section II. A and Figure 1). 
The statistical model of erosion rate in its present form supposes that, 
at the start of a test run, all of the a-cells in the original surface 
are equivalent. In order to treat high-strength materials successfully, 
it will be necessary to restrict the number of cells that are subject to 
crack formation as the result of a drop impact. An attempt was made to 
apply the present form of the statistical model to aged Udimet 700 alloy; 
the obstacles that were encountered are described below. 

Because of the exceptionally long test times that were required to 
erode Udimet 700 alloy, it was necessary to collect all of the volume- 
loss data on one test specimen. This introduced reset errors in the 
impingement area, A, as a consequence of removing and replacing the test 
specimen after each increment of test in order to determine the weight 
loss. The reset error can only be an increase in the impingement area, A. 
In calculating the probability of a hit, p, the first number of the count- 
ing rule, N 0 , and the area of a microcell, the impingement area was arbi- 
trarily assumed to have increased from 0.0235606 cm^ to 0.0352991 cm^; 
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Figure 14. Theoretical Rate-of- Volume-Loss and Ac cumulated- Volume-Loss Curves for Zinc 
Calculated with Fragments of Equal Volume. 


Accumulated Volume Loss 




this is essentially a 50 percent increase. Such an increase in the im- 
pingement area has only a small effect on the locus of the theoretical 
curves. This is shown in Section V. 

The experimental volume-loss data for Udimet 700 aged do not show 
the marked step structure that is apparent in the experimental volume- 
loss data for iron, nickel, and zinc. Nevertheless, some structure 
exists; this can be seen by inspection of the aged Udimet 700 volume- 
loss data shown in Figure 15. From Figure 15 it can be seen that there 
are sections of the plot that can be fit by straight lines and that these 
lines gradually change in slope. 

From inspection of the experimental data plotted in Figure 15, the 
first sharp rise in accumulated volume loss was taken to occur at 
n 0 = 2.2 x 10 impacts. On the basis of this value of n Q , the proba- 
bility of a hit against the typical cell was found to be p = 0.0000011236 
and the first number of the counting rule, N Q , was found to be 2.2472. 

The fact that the first number of the counting rule is less than 
three presents an anomaly because a minimum of three impacts has been 
accepted as necessary to circumscribe a piece of the solid which then 
breaks away as an eroded fragment. A number less than three for N q can 
only be justified by accepting it as a statistical average and by assuming 
that residual crack ends contribute to the removal of new fragments. 

The diameter of an eroded aged Udimet 700 fragment was found to be 
2.2472 microns. This is about 1/70 of the average grain size of the aged 
Udimet 700 metal used for the test specimen; the average grain size was 
found to be 160 microns with use of the Heyn intercept method. 

The rate of volume loss and the accumulated volume loss found for 
aged Udimet 700 with use of the computer program for Equation 7 are 
plotted against number of impacts sustained in Figure 16 along with the 
experimental data points. From this plot it can be seen that the theo- 
retical rate curve is quite different from that found for the other metals 
and that the accumulated-volume-loss curve is in agreement with the ex- 
perimental data only over the range from 2 x 10^ to 5 x 10^ impacts. 

Better agreement is found when the change in volume of eroded fragments 
is taken into account (see Section IV) but the calculated values of the 
probability, p, and of the first number of the counting rule, N 0 , for 
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Figure 15. Zones in the Experimental Volume-Loss Data for Aged Udimet 700. 



Rate of Volume Loss, ec/impact 



Accumulated Volume Loss, cc 







Udimet 700 aged suggest that this metal is at or beyond the limits of 
valid use of the cellular model of erosion rate in its present state of 
development. 

G. Features of Rate Curves 

As far as drop-impact erosion is concerned, solids may be classified 
as Class A materials or Class B materials (see Section II. A) depending on 
whether or not they fail as the result of the radial tensile stresses that 
are imposed by each drop impact that occurs under an arbitrary set of 
erosive conditions (relative impact velocity, drop size, and density of 
drop liquid). The first approach to a cellular model of erosion rate, 
which has been described in the preceding sections of this report, is 
directed to Class A solids, that is, to solids that fail under the tensile 
stresses imposed by each drop impact. 

Because of the geometry of fracture (see Section II. B), the erosion 
of a Class A solid will be characterized by the removal of layers of 
material. The erosion-rate feature associated with the removal of a 
layer of material from a solid under drop-impact test is a pulse in which 
the erosion rate increases to a maximum value and then decays to zero. 

The calculated rate curves for iron, nickel, and zinc found with use of 
Equation 7 consist of a series of such pulses; these pulses correspond 
to the removal of layers of material from the test specimens. The ex- 
perimental volume-loss data are in agreement with the calculated rate 
curves because iron, nickel, and zinc are Class A solids with respect to 
the erosive conditions that were imposed in collecting the experimental 
volume-loss data. 

The statistical model has not yet been developed to include Class B 
solids which do not fail under the radial tensile stresses imposed by 
each drop-impact; Class B solids fail only at weak spots after a large 
number of impacts has been sustained (see Section II. A). It seems reason- 
able to expect that the calculated erosion-rate curves of Class B solids 
will not consist of a series of pulses. This expectation is based on the 
reasoning that Class B solids will only fail after a large number of 
impacts has been sustained and that the probability distribution for a 
hit in a a-cell will spread out as the number of impacts delivered is 
increased (see Section II.E.l.c). It is inadvisable to say more about 
the form that the calculated rate curves of Class B solids may have until 
the statistical model has been developed to include materials of this class. 
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An attempt to make a partial test of the equations with experimental 
volume-loss data for Udimet 700 alloy (see F in this Section) has provided 
evidence which suggests that this alloy is outside the limits of applica- 
bility of the statistical model in its present state of development. The 
very high strength of Udimet 700 alloy in comparison with the strengths 
of iron, nickel, and zinc strongly suggests that Udimet 700 alloy is a 
Class B solid with respect to the erosive conditions that were used. 

It is important to recognize that whether a solid is in Class A or 
in Class B depends both on its characteristic strength properties and/or 
ability to accept energy elastically and on the severity of the erosive 
conditions that are used in testing it. As in all attempts to separate 
classes, there will inevitably be an intermediate or transition class. 

The features of the erosion-rate curve of a given solid will be different 
depending upon whether it is Class A, Transition Class, or Class B with 
respect to the severity of the erosive conditions to which it is subjected. 
The erosion-rate curve found for tantalum (see D in this Section) suggests 
that tantalum may really be in the Transition Class with respect to the 
erosive conditions that were used because the pulses in the rate of 
volume loss do not drop to zero as they do in the case of iron, nickel, 
and zinc . 

H. Correlation Between Quantities N Q and p and Physical Properties of 
Metals 

From the time of the earliest investigations of drop-impact and 

cavitation erosion it has been found that erosion resistance correlates 

in general with the tensile strength and hardness of the metal tested. 

This suggests that correlation should be sought between quantities that 

appear in Equation 7 and the tensile strength and hardness of the selected 

metals. The probability of a hit, p, and the first number of the counting 

rule, N Q , which were found for the selected metals with recourse to the 

experimental volume- loss data are listed in Table 3 along with the 0.2 

percent offset yield strength, the ultimate tensile strength, and the 

Knoop microhardness number of the particular heat and heat-treatment 

(14) 

state of each selected metal. 

It can be seen by inspection of the tabulated data that inverse 
trends exist between the value of the probability, p, and both the 0.2 
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TABLE 3 


A. QUANTITIES THAT APPEAR IN EQ (5) AND PHYSICAL PROPERTIES OF THE SELECTED METALS 


Metal 

Probability, p 

N 

o 

0.2% Offset 
Yield Strength 

Ultimate 

Tensile Strength 

Knoop 

Microhardness No. 


— 

impacts 

dyn/cm^ 

2 

dyn/cm 

— 

Zinc 

0.0300018 

300 

a 

4.5 

X 

00 

o 

H 

a 

1.07 

X 

io 9 

3 32 

Nickel 

0.00333313 

100 

a 

6.3 

X 

o 

00 

a 

3.54 

X 

io 9 

a 57 

Iron 

0.00043135 

32 

a 

14.8 

X 

00 

o 

iH 

a 

2.76 

X 

io 9 

3 72 

Tantalum 

0.000208321 

25 

a 

15.7 

X 

00 

O 

rH 

a 

2.40 

X 

io 9 

3 104 

U-700 Ag. 

0.0000011236 

2.2 

a 

89.7 

X 

00 

o 

a 

13.8 

X 

io 9 

3 375 


B. OTHER EROSION-RATE QUANTITIES AND STEEL-SPHERE DENTING VELOCITY 


Metal 

Maximum Rate 

Impacts to 

Eroded 

Average 

Steel Sphere 


Volume Loss 

Reach Maximum 

Fragment 

Grain 

Denting 



Rate Vol. loss 

Volume 

Size 

Velocity 


3 

cm /impact 

impacts 

3 

cm 

micron 

cm/sec 

Zinc 

38.8 x 10 -8 

10,000 

366 x 10 8 

300 

3 15.8 

Nickel 

-8 

9.6 x 10 

30,000 

17.6 x 10~ 8 

53 

a 11.2 

Iron 

-8 

1.8 x 10 

91,000 

0.58 x 10~ 8 

63 

a 

35.3 

Tantalum 

-8 

0.86 x 10 

120,000 

0.20 x 10~ 8 

126 

a 38.6 

U-700 Ag. 

0.061 x 10~ 8 

2,100,000 

0.035 x 10 -8 

160 

3 64.2 


Values measured for the specific heat and heat-treatment state of the selected 
metals. See Reference 14. 

Values found for average grain size with use of the Heyn intercept method. 
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percent offset yield strength and Knoop microhardness number; a corres- 
ponding trend does not exist through the values of ultimate tensile 
strength. A plot of negative log p against 0.2 percent offset yield 
strength is given in Figure 17. In drawing the line shown -in this 
figure, correlation was sought principally for iron, nickel, and zinc, 
which are Class A metals with respect to the erosive conditions that 
were used. From Figure 17 it can be seen that tantalum, which may be 
a Transition-Class metal with respect to the erosive conditions used, 
also fits on the line. Udimet 700 alloy, which is a Class B metal with 
respect to the erosive conditions used, does not fit on the line. 

An inverse trend exists between the first number of the counting 
rule, N 0 , and both the 0.2 percent offset yield strength and Knoop 
microhardness number. Plots of log N Q against 0.2 percent offset yield 
strength and Knoop microhardness number are shown in Figures 18 and 19, 
respectively. Correlation was sought principally for iron, nickel, and 
zinc; the difference between tantalum and the purely Class A metals 
iron, nickel, and zinc shows up in Figure 19. Udimet 700 alloy, which 
is a Class B metal with respect to the erosive conditions used, shows 
no correlation with the other metals in either plot. 

From Table 3B it can be seen that an inverse trend exists between 
both the maximum rate of volume loss and the eroded fragment volume and 
the yield strength of the metals. On the other hand, the number of 
impacts required to reach the maximum rate of volume loss varies directly 
with the yield strength of the metals. 

I. Experimental Confirmation of the Number of Fragments in an Eroded 
Layer 

An effort was made to obtain experimental confirmation of the intui- 
tive argument that an eroded layer of material contains nn 1 fragments by 
calculating the volume of material lost when the test specimen was 
punctured by erosion or was near the point of puncture. The numbers of 
impacts for which this volume loss was calculated for each metal are 
given in column two of Table 4. The computer print-out of layers re- 
moved, which is designated as Z, is given for each of the metals in 
column three. Values of the product nn' for the metals are given in 
column four and the numbers of fragments removed, taken to be the 
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Figure 17. Correlation Between the Probability of a Hit, p, and the 0.2 Per Cent Offset Yield Strength. 








TABLE 4 

CALCULATED VOLUME OF METAL REMOVED FROM A TEST SPECIMEN BY LIQUID DROP IMPINGEMENT 


Metal 

Impacts , 
n 

Layers 
Removed , Z 

on’ 

Fragments 
Removed, ' Z 

Volume per 
Fragment , cc 

Volume of Metal 
Removed, cc 

Impacts Actually 
Sustained by the 
Test Specimen 


5 




, „-9 


5 

Tantalum 

15 x 10 

4.76766 

261209. 

1245355. 

1.96997 x 10 

0.002453 

15.5 x 10 








Hole Formed 


5 




-7 


5 

Nickel 

4.4 x 10 

3.19425 

4081.41 

13037.1 

1.75596 x 10 

0.002289 

4.36 x 10 


5 




-9 


5 

Iron 

10 x 10 

3.28918 

124532 

409609. 

5.79255 x 10 

0.002373 

9.69 x 10 


5 




-6 



Zinc 

0.4 x 10 

3.50484 

151.137 

529.711 

3.65924 x 10 

0.001938 

0.415 x 10 








Hole Formed 



product nn'Z, are given in column five. The initial volume of an eroded 
fragment for each of the metals as found by trial (see Section B through 
F) is given in column six and the volume of metal removed, taken to be 
the product of the number of fragments removed multiplied by the volume 
of a fragment, is given in column seven. The actual number of impacts 
sustained by a test specimen of each of the metals is given in column 
eight; for two of the metals (tantalum and zinc) erosion progressed 
through the thickness of the specimen. 

The thickness of the test specimens, as determined by a measurement 

on one specimen only, was 0.33 cm. The area over which impingement 

occurred was calculated to be 0.02356 cm^ on the basis of the observation 

that the impinging drops struck over an area the radius of which was 

roughly equal to the drop diameter. The volume of the cylinder of metal 

that would have been lost had this area been moved through the thickness 

3 

of the test specimen is 0.00778 cm . It can be seen that the calculated 

volume of metal removed given in column seven of Table 4 is equal to 

3 

about 30 percent of the quantity 0.00778 cm , which is the total possible 
volume of metal that could have been removed. 

On one hand, the total possible volume of metal that could be 

3 

removed, assessed as 0.00778 cm on the basis of removing a cylinder of 
metal with a base area equal to the impingement area and a height equal 
to the thickness of the test specimen, is too high because erosion 
craters have rounded bottoms (see Figure 20) . On the other hand, the 
calculated volume of metal removed is too low because the volume of 
eroded fragments increases as tests progress. It is considered that the 
extent of agreement found is sufficiently good to confirm the use of 
product rm' as a multiplying factor to obtain fragments removed from 
layers removed. 
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A 


Hole through a plate formed by removal 
of a cylinder of the plate metal 



B 


Hole through a plate formed by the 
growth of an erosion crater 


Figure 20. Evidence to Show that the Volume of Metal Removed by the Growth 

of an Erosion Crater is Less Than That Removed by Moving the Area 
of the Mouth of the Crater Through the Thickness of the Test 
Specimen. 
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IV. ACCUMULATED VOL U ME LOSS FOUND WITH FRAGMENTS OF CHANGING SIZE 


The concept of layers of fragments of constant volume given as 
integral multiples of nn' in Equation 7, can be expanded to associate 
the first layer of eroded fragments with a volume v^, two layers of 
eroded fragments with a volume v^ + v 2> nnd three layers of eroded 
fragments with a volume v^ + v 2 + V 3 * v: *- ew °f these changes. Equation 
7 becomes 


L(x,n) = v 1 


N + N. - 1 
o 1 

xz 


P(x,n) + (v 1 + v 2 ) 


N + N. + N„ 
O 1 z 

YZ 


- l 


P(x,n) 


x = N 

o 


x = N + N. 
o 1 


+ (v x + v 2 + v 3 ) 


N + N n + N„ + N„ - 1 
o 1 l 3 


P(x,n) + ... (18) 


x = N + PL + N_ 
o 1 l 



where L(x,n) is now the accumulated volume loss found directly for a 
given number of impacts, n. For the case that v^ = v 2 = v^. Equation 18 
is as well founded as Equation 7. It has the advantage that accumulated 
volume loss is printed out directly by the computer program; it has the 
disadvantage that erosion rate is not obtained as an intermediate step 
in the calculation of accumulated volume loss. 

For the case that v^ # v 2 t v^. Equation 18 is not as well founded 
as Equation 7. If the equality between cell size and fragment size is 
maintained and both the cell size and fragment size are allowed to in- 
crease (or decrease) during a calculation, a statistical problem is 
encountered; a corresponding change in the probability of a hit must be 
made. On the other hand, if the cell size and, consequently, the proba- 
bility of a hit are maintained constant throughout a calculation but the 
volume of a fragment or of a layer of fragments is allowed to change, 
then a physical problem is encountered; the radius of the circle of 
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fracture is fixed because it is determined by the drop size and it is 
not clear how a change in average fragment size can be accounted for 
(see Figure 4) if the fragment dimensions are increased (or decreased) 
equally. If, however, it is assumed that the fragments increase (or 
decrease) in volume because the depth at which the cracks intersect to 
release a fragment becomes greater (or less), no problem is encountered. 
With recourse to this argument, the statistical model and the computer 
program in the state of development in which they now exist can be used 
to make calculations in which the volume of a fragment changes. 

Equation 18 was used to calculate accumulated-volume-loss curves of 
the selected metals for the case that the volume of a layer of fragments 
changes. For each metal, the same probability, p, of a hit against the 
typical cell and the same counting rule, N^, were used as were used in 
obtaining the calculated curves of Section III. The value of v^ was 
assessed directly from plots of experimental volume loss of the metals 
plotted against number of impacts sustained. It was taken to be the 
volume loss reached at the first knee in these plots. Values of the 
succeeding v^ were also assessed approximately from inspection of the 
experimental volume-loss plots; they were expressed in terms of v^ for 
each of the metals. The direct printout of the accumulated volume loss 
for each of the five selected metals was plotted against the numbers of 
impacts sustained. The plots obtained are shown in Figures 21, 22, 23, 

24 and 25. Comparison of Figures 21-25 with Figures 11-14 and 16 shows 
that in almost every case the correction for increasing fragment size 
is all that is needed to establish good agreement between the theoretical 
curves and the experimental data. 

The volume of the first layer, v^, and the nature of the percentage 
layer-volume change for the selected metals is given in Table 5. It can 
be seen that in the case of aged Udimet 700, and eventually in the case 
of tantalum, the change in layer volume is a decrease . Because Udimet 
700 is a Class B metal and because tantalum may be Transition-Class metal 
with respect to the erosive conditions that were used, this observation 
suggests that Class B metals may be characterized by a decrease in eroded 
fragment size. The 0.2 percent offset yield strength and the ultimate 
tensile strength of each metal is also given in Table 5. There appears 
to be no correlation between the 0.2 percent offset yield strength and 
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Volume Loss , cc 


o 


Experimental Data Points □ 

Input Data for Calculated Curve: 
Probability, p = 0.000208321 
Counting Rule, N: N = 24.9985 impacts 

= 49.9970 impacts 
N£ = 74.9955 impacts 
Volume of First Layer: = 0.0005 cc 

Vj = 2 = 0.0010 cc 

v„ = 0.0004 cc 


Number of Impacts 


Figure 23. Calculated Accumulated-Volume-Loss Curve for Tantalum with 
Fragments of Increasing Volume. 
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TABLE 5 


LAYER VOLUME AND PERCENTAGE CHANGE IN LAYER VOLUME 


o 


Metal 

Layer Volume, v 

Percentage Change 
in Layer Volume 

0.2% Offset 
Yield Strength 

Ultimate 

Tensile Strength 


3 

cm 


psi 

psi 

Udiraet 700 

0.0020 

62.5% Decrease 

130,000 

201,000 

Nickel 

0.00072 

10% Increase 

9,000 

51,000 

Iron 

0.00072 

20% Increase 

21,000 

43,000 

Tantalum 

0.0005 

100% Increase then 20% 
Decrease 

23,000 

35,000 

Zinc 

0.0006 

50%> Increase 

7,000 

15,000 










the layer volume, v^, or the percentage change in layer volume. However, 
the ultimate tensile strength appears to correlate directly with the 
layer volume, v^, and inversely with the percentage change in layer 
volume (see Figures 26 and 27) . 

In Figures 26 and 27 the volume of the first layer, v^, and the 
percentage change in layer volume in terms of v^ are plotted, respectively, 
against the ultimate tensile strength. Correlation was sought principally 
for iron, nickel, and zinc, which are Class A metals with respect to the 
erosive conditions that were used. The slight divergence of tantalum and 
the wide divergence of Udimet 700 alloy are evident in these plots. These 
divergences reflect the fact that Udimet 700 is a Class B metal and that 
tantalum may be a Transition-Class metal with respect to the erosive 
conditions that were used. 
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Percent Change in Layer Volume 






V. EFFECT OF A CHANGE IN SIZE OF THE IMPINGEMENT AREA. A 


To what extent a change in size of impingement area. A, may affect 
the locus of the accumulated-volume-loss curve is an important consider- 
ation. In collecting drop-impact volume- loss data, area A was assessed 
visually by watching drops impinge and comparing the diameter of the 
area over which impingement occurred with the diameter of a drop. The 
value of area A assessed in this way is not known with a high degree of 
reliability. It is, therefore, important to know to what extent agree- 
ment between the theoretical accumulated-volume-loss curve and the ex- 
perimental data may be affected by a change in the value of area A. An 
additional reason for determining the importance of a change in the 
value of area A on the theoretical accumulated-volume-loss curve is that 
in collection of cavitation erosion data with use of a magnetos trictive 
device the eroded area tends to increase in size during a test run. 

The accumulated-volume-loss curve for iron was calculated with an 
approximately 50 percent increase and 50 percent decrease in the size of 
area A. The calculation was carried out with use of Equation 18 of 
Section IV. The calculated curves are plotted along with the experimental 
accumulated-volume-loss data for iron in Figure 28. The curve for an 
impingement area A equal to 0.0235606 in Figure 28 is the same curve as 
that given for iron in Figure 21. 

Two conclusions can be drawn on the basis of the curves shown in 
Figure 28. The first conclusion is that changing the size of the im- 
pingement area by a factor of three is not important as far as the fit 
of calculated curves with experimental data points is concerned for 
metals with an erosion resistance comparable to that of iron. There is 
reason to suppose that this effect may be more important for metals that 
are highly erosion resistant. The second conclusion is that the step- 
structure of the calculated accumulated-volume-loss curve is sharper and 
better delineated for large than for small values of the probability, p, 
and of the counting rule, N^. (See data inset in Figure 28). On the 
basis of the relation given by Equation 6, an increase in size of area A 
should have the effect of smoothing out the step structure of a calculated 
curve. 
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Accumulated Vol 



Impingement 
Area, A, in 
sq. cm. 

0.0353409 

0.0235606 

0.0117803 

__J 

6 


Counting Rule, 


15.9815 31.9630 

32 64 

47.9431 95.8862 


N 


63.9260 — 
128 _ 
191.7724 


7 8 9 10 x 10 


5 


Impacts 

Over Which Impingement Is Occurring. 




VI. CONCLUDING NOTE 


The work on a statistical model of erosion rate described in this 
report is by no means finished. Two quantities, namely, n Q , the number 
of impacts that must occur against area A before the first sharp rise in 
volume loss occurs, and v, the volume of an eroded fragment, can at 
present only be assessed by recourse to experimental data for volume loss 
with elapsed time or number of impacts sustained. Although the relation 
between impact energy and velocity and the volume of an eroded fragment 
has been established for two shapes of fragments, no work has as yet been 
done to determine the dependence of the volume of an eroded fragment on 
the strength of the solid material against which impingement occurs. 

In addition to completing the statistical model for Class A brittle 
materials which fail under each drop impact, the model needs to be de- 
veloped further to include the Class B brittle materials whose ability 
to accept energy is large in comparison with the impact energy delivered 
by each drop blow. Finally, the model should also be extended to include 
materials that are permanently ductile or plastic and even rubbery materials 
such as neoprene. What has been presented in this report is really only 
a beginning. 


76 



VII. REFERENCES 


1. Naude, C. F. and Ellis, A. T., Trans. A.S.M.E. 83, 648, 1961. 

2. Thiruvengadam, A. Preiser, H. S., Journal of Ship Research J3, 39, 

1964. 

3. Hammitt, F. G.; Robinson, M. J.; Seibert, C. A.; and Aydinmakine, 

F. A.; University of Michigan Report 03424-14-T (NASA Grant No. 

NsG— 39-60) , 1964. 

4. Heymann, F. J., Symposium on Erosion and Cavitation, A.S.T.M* 

Special Technical Publication No. 408, 1967. 

5. Mok, C. H., Mechanics Section TM #068-11, General Electric Company 
Space Sciences Laboratory, Philadelphia, Pennsylvania, 1968. 

6. Hammitt, F. G. ; Huang, Y. C. ; Kling, C. L.; Mitchel, T. M. ; and 
Solomon, L. P. ; University of Michigan College of Engineering Report 
No. 02643-2-PR, 1969. 

7. Hoff, G.; Langbein, G. ; and Rieger, H. ; Symposium on Erosion and 
Cavitation, A.S.T.M. Special Technical Publication No. 408, 1967. 

8. Engel, 0. G., Discussion of a Paper given by F. J. Heymann, A.S.T.M. 
Special Technical Publication No. 408, 1967. 

9. Gould, G. C., Proceedings of the A.S.T.M. Symposium on Character- 
ization and Determination of Erosion Resistance, Atlantic City, N.J., 
1969. 

10. Schwarz, M. v. and Mantel, W., Zeit. des Vereines Deutscher Ingenieure 
80, 863, 1936. 

11. Baker, D. W. C.; Jolliffe, K. H. ; and Pearson, D. ; Phil. Trans. Roy. 

Soc. 260A , 193, 1966. 

12. Engel, 0. G., National Bureau of Standards Jour, of Research 64A, 

499, 1960. 

13. Engel, 0. G., Wright Air Development Center, TR 53-192 Part III, 1953. 

14. Engel, 0. G., General Electric Company, Report GESP-253 (Contract 
NASW 1481), 1969. 

15. Engel, 0. G. , Wright Air Development Center, TR 53-192 Part VII, 1957. 

16. See, for example, Engel, 0. G. , National Bureau of Standards Jour, of 
Research, _54, 51 (1955) and Brunton, J. H., A.S.T.M. Special Techni- 
cal Publication No. 307, 1961. 

17. Auleytner , J. , X-Ray Methods in the Study of Defects in Single Crystals , 
Pergamon Press, New York, 1967. 


77 


I 


REFERENCES (Cont'd) 


18. Barrett, C. S., Structure of Metals . McGraw-Hill Book Co., Inc., 

New York, 1943. 

19. Friedel, J. , Dislocations . Addison-Wesley Publishing Co., Inc. 
Reading, Massachusetts, 1964. 

20. Bragg, W . L., Proc. Phys. Soc. 52 , 54, 1940. 

21. Burgers, J. M. , Proc. Phys. Soc. J52, 23, 1940. 

22. Brandenberger , E. and DeHaller, P., Schweiz. Arch. 10, 331 and 379, 
1944. 

23. Hoff, G. and Langbein, G., Proc. Second Meersburg Conference on Rain 
Erosion and Allied Phenomena, Vol. 2, Ministry of Technology, Royal 
Aircraft Establishment, Farnborough, Hants, England, 1967. 

24. DeHaller, P. , Schweizerische Bauzeitung 101 . 243 and 260, 1933. 

25. See, for example, Kern, W. F. and Bland, J. R. , Solid Mensuration , 
2nd Ed., John Wiley and Sons, Inc., New York, 1938. See page 51. 

26. Griffith, A. A., Phil. Trans. Roy. Soc. 221A. 163, 1920-21. 

27. Adamson, A. W. , Physical Chemistry of Surfaces , Interscience 
Publishers, Inc., New York, 1960. 

28. Ripken, J. F., A.S.T.M. Special Technical Publication No. 408, 1967. 


78 



VIII. LIST OF PRINCIPAL SYMBOLS 


A 

a 

r 

n 

n* 


d 

N. 

1 

N 

o 

V 

M 

e 

R 

J 

v 

n 

P 


P(x,n) 

F(x,n) 


Q 

L 

n 

o 


v 

c 



area over which impingement is occurring 

area of a single cell in the array of cells covering area A 
radius of the drops that impinge 
A/a = number of cells in area A 

number of cells that develop cracks when a single drop impact 
occurs 

diameter of a cell and also distance across an eroded fragment 
numbers of the counting rule 

first number of the counting rule or number of impacts required 
to remove a layer of fragments from the original surface 

impact velocity 

drop mass 

angle of attack measured from the surface of the test specimen 
rate of volume loss 

average rate of ejection of fragments from area A 

average volume of an eroded fragment 

some arbitrary number of impacts against area A 

a / A = probability of a hit against any one of the cells in 
area A 

probability of x impacts on cell area, a> when n impacts have 
occurred on area A 

number of fragments ejected from area A after n impacts have 
occurred on area A 

incremental volume loss 

accumulated volume loss 

r|N = number of impacts that must occur against area A before 
the first sharp rise in volume loss occurs 

minimal volume of an eroded fragment 

volume of a prismatic eroded fragment 

volume of a tetrahedral eroded fragment 


79 


LIST OF PRINCIPAL SYMBOLS (Cont'd) 


E = the part of the impact energy invested in crack formation 

E c = amount of energy needed to form a crack of smallest size or 
the minimal value of E 

V = impact velocity for which the energy invested in crack formation 

is E c 

A c = area of an eroded platelet of smallest possible size 
t CT = thickness of an eroded platelet 

T = surface tension 

Y = Young’s modulus of elasticity 

v = Poisson's ratio 

c = half length of a Griffith crack 

Sy =0.2 percent offset yield strength 
S = strength of a solid 
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APPENDIX A 


POSSIBILITY OF IMPACT IONIZATION OF WATER 


The heat of neutralization of water at infinite dilution is 13.4 
Kcals/mole, that is, 

H + (liq) + OH~(liq). ^ ' — H 9 0(liq) AH = 13.4 Kcals/mole. 

The volume of a 0.2-cm drop is (4/3) (3.1416) (0.001) or 4.2 x 10 - ^ cm^. 
For water, this drop volume is equivalent to 4.2 x 10"^ grams. Because 
one mole of water is 18 grams, the heat of neutralization per drop is 
found to be (4.2 x 10~^/18) [l3.4]or 3.1 x 10 - 3 Kcals. 

If the relative collision velocity is 880 ft/sec (26,822 cm/sec), 
the kinetic energy of the impact is (1/2) (4.2 x 10“ 3) (2.682 x 10^)2 or 
3.6 x 10 - ^ Kcals. This is about 0.01 of the heat of neutralization per 
drop, which was found to be 3.1 x 10~3 Kcals. In terms of this result 
it is reasonable to expect that some ionization will occur. At relative 
collision velocities higher than 880 ft/sec, ionization will become more 
probable. 
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APPENDIX B 


COLLECTION OF DROP-IMPACT-EROSION DATA 
by John A. Almo 

I. EXPERIMENTAL APPARATUS 

The basic facility consists of a 24 inch diameter rotor with a test 
specimen attached at the periphery, rotating in a horizontal plane. A 
vertical stream of water drops is generated and made to fall into the 
path of the rotating test specimen. Photoelectric cells triggered by a 
slotted disk connected to the rotor provide signals for synchronized 
drop generation and for stroboscopic observation of the rotor and falling 
drop stream. Additional electronic apparatus is used to provide constant- 
speed operation at any target velocity from 400 ft/sec to 1250 ft/sec, to 
count the number of drop impacts and to control the relative position of 
the drop stream and the rotating specimen. (Figures 1 and 2 show the 
rotor assembly with a 1/4-inch diameter test specimen mounted in the 
rotor tip.) 

The rotor is enclosed in a protective steel chamber fitted with 
thick plastic viewing windows. The chamber is evacuated during oper- 
ation to reduce power requirements and eliminate windage disturbance of 
the drop stream. (Figure 3 shows the test rig with the protective chamber 
closed. A view of the test rig with the chamber opened to permit in- 
spection or change of the test specimen is shown in Figure 4.) 

The rotor system, essentially patterned after those used in ultra- 
centrifuges, provides extremely stable high-velocity operation. The 
test apparatus can be operated continuously for up to four hours during 
which time water drops strike the target with dispersion no greater than 
two drop diameters. 

II. APPARATUS CHARACTERISTICS AND TEST PROCEDURES 

The following are descriptions of several apparatus characteristics 
and test procedures which have a bearing on the reliability and inter- 
pretation of the test data obtained. 
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A. Water Drops 

Degassed, distilled water was used throughout the test program. 

The water drops (generated by a vibrating capillary tube) experienced 
oscillatory shape changes immediately after formation. These oscil- 
lations were rapidly damped out in the vertical fall from the drop 
generator to the impact target, however, and drop shape at the time of 
impact was observed to be essentially spherical. 

Drop volume was derived from measurements of flow rate to the drop 
generator and of drop generation rate. This method of measurement gave 
a value of drop volume accurate to ±1%. Drop diameter was then calculated 
assuming a spherical drop of the measured volume. 

B. Target Velocity 

Target velocity was calculated from measurements of the rotor 
diameter and rotation rate. The rotation rate was monitored continuously 
throughout the test period and the impact process was interrupted if the 
target velocity drifted outside the limits ±1 ft/sec. 

C. Drop Alignment and Impact Count 

Drop generation was initiated after stabilization of the rotor 
speed. Alignment of the drop stream and the target specimen was in- 
directly accomplished with the drop stream descending clear of the 
rotor tip. Once good alignment was obtained externally, the drop 
stream was electrostatically deflected into the path of the target. 
Simultaneously, an electronic counter was connected to the deflection 
circuit to measure the total number of deflected, and thus impacted, 
drops. The synchronization of the drop generator, deflection circuit 
and rotor was such that an impact with a deflected drop occurred for 
every rotor revolution. 

The operator continuously observed the target specimen through a 
telescope and allowed the count to proceed unless he observed excessive 
drift in the impact location. If this occurred, he could interrupt the 
drop deflection and realign the drop impact point. Once realigned, the 
deflection was reestablished and the impact count continued. 

For run times exceeding five minutes the elapsed time on the target, 
rather then the number of impacts, was measured. For convenience, all 
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data are listed in elapsed time, but may easily be converted to number 
of impacts using the known impact rate. 

D. Impinged Area 

At the 1000 ft/sec operating speed used in these tests, windage 
effects from the remaining air in the test chamber produced a random 
impingement location circumscribed by a circle of diameter equal to two 
drop diameters. All tests, with two exceptions to be noted later, were 
carried out under these dispersion conditions. 

E. Drop Splash 

Observations of the conical spray pattern produced at impact suggest 
that the initial lateral liquid flow following impact does not occur in 
a purely radial direction, but possesses an additional component normal 
to the specimen face. No quantitative measurements of this flow were 
made. 

F. Ambient Test Pressure 

Throughout the test program the ambient pressure was maintained at 
0.02 atmospheres. No evidence of drop instability (indicative of boiling) 
was observed. At pressures greater than 0.02 atmospheres, the windage 
disturbance of the drop streams increased, resulting in an increased 
dispersion of the impact location. At pressures lower than 0.02 atmo- 
pheres, drop stability was adversely affected. 

G. Ambient Test Temperature 

The temperature on all walls of the test chamber was kept at 10°F 
to condense water vapor formed by the impact process. The drop generator 
expelled drops at essentially room temperature, but the drop and target 
temperatures at the time of impact are unknown. Tests conducted at 
speeds below 500 ft/sec in an earlier program showed an ice formation 
around the target area. At 1000 ft/sec no ice formation was evident, 
but radial forces were probably too great to allow ice formation. 

H. Specimen Handling 

All targets were initially inspected for obvious surface chips, 
dents, or other mechanical irregularities and individually washed in 
methyl ethyl ketone to remove machine oil and grease. After the MEK 
wash, the targets were handled only with forceps or surgical gloves to 
avoid further oil contamination. The targets were then thoroughly washed 
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in distilled water, dried for one hour in a hot air stream, and stored 
in a desiccator jar in individual specimen envelopes. 

After 24 hours in the desiccator, the targets were weighed on a 
laboratory balance (accuracy ± 0.05 milligram) and returned to storage 
in the desiccator. As a precaution, the targets were weighed again just 
prior to testing and any discrepancies with the first measurement were 
resolved by additional weighings before the targets were mounted on the 
rotor arm for testing. Immediately after testing the targets were 
gently washed in distilled water, dried in the hot air stream, and re- 
turned to the desiccator. The final weights were then recorded after 
24 hours in the desiccator. Again as a precaution each target was 
weighed at least twice. 

I. Test Methods 

The weigh t-loss-versus-exposure- time data for zinc, iron, tantalum, 
and nickel were obtained by a single- test— point-per-sample method. Each 
data point represents a different sample of the material tested for the 
indicated time. This method of testing was chosen so that individual 
specimens representing the various exposure times could be retained for 
physical examination. The Udimet materials were tested in the usual 
cumulative-weight-loss-on-a-single-sample manner for reasons to be ex- 
plained later. 


m . RESULTS 

Six materials were tested: zinc, iron, tantalum, nickel, Udimet 700 

(aged), and Udimet 700 (solutioned) . A discussion of the tabulations of 
weight loss versus exposure time for each material is presented below. 

A. Zinc (see Table 1) 

Little or no "incubation" time was observed for the zinc samples 
under the conditions of testing. The shortest exposure time taken, 

Sample 2, showed a very distinct crater after 15.2 sec of exposure. 

Distinct ridges around the impact crater are probably indicative of 
plastic deformation. These were observed on Samples 2, 4, 5, 10, 12, A, 
and C. Samples 14 and 18 are not reliable, since a hole was punched 
through the target during the test run and the time at which the 
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punch- through occurred could not be observed. During the test run on 
Sample 13 the impact location drifted vertically approximately 5 drop 
diameters; Sample 13 is therefore also unreliable. 

Sample 15 was hit at two separate locations. At one a small crater 
formed after 11.9 sec and at the other a much larger crater formed after 
180 sec of impact time. The final weight figure shown for this sample 
is the sum of the losses for the two craters; the weight loss for the 
11.9 sec exposure crater alone was 0.2 milligram. 

B. Iron (See Table 2) 

The distinct ridge around the crater evident in many of the zinc 
samples was present only on Sample 1 of the iron series. The additional 
column labeled "initial damage" lists the exposure time at which the 
operator first observed damage on the test specimen face. Althouglji the 
test procedure from sample to sample was unchanged, a wide variation in 
initial damage time was observed. 

C. Nickel (See Table 3) 

Samples 1, 8, and 12 of the nickel series had distinct ridges 
surrounding the impact crater. The nickel tests were conducted early 
in the program before it became apparent that initial damage time might 
vary significantly from sample to sample. For this reason only four 
values of initial damage time were recorded. 

Unlike the other materials tested, most of the nickel samples 
showed clear evidence of some form of secondary impact damage. With 
the exception of Samples 1 and 2, which had the lowest exposure times, 
the nickel targets were characterized by the appearance of small dents 
distributed over the face of the target after an impact erosion crater 
had formed. The dents did not have the same shape as the primary crater 
formed at the impact point, but rather appeared to be small material de- 
formations with little or no actual material loss. Possibly these 
indentations were formed by the impact of the target with small material 
fragments removed from the drop impact crater. Since only the nickel 
specimens exhibited the secondary impact damage, the secondary impact 
of the target with the liquid spray cone formed at the drop impact point 
(which occurred for all materials tested) is probably not the damage 
mechanism. 
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D. Tantalum (See Table 4) 

Samples 1, 3, and 8 had distinct ridges surrounding the impact 
crater. Sample 2 was inadvertently run at a pressure of 0.04 atmospheres, 
resulting in an impact area 3 drop diameters wide. Sample 9 was punched 
through sometime during the test run and is therefore unreliable. The 
initial damage times observed had large variations, again with no vari- 
ations in test procedure from sample to sample. 

E. Udimet 700 (See Tables 5 and 6) 

The U-700 materials proved to be very resistant to damage and were 
therefore tested for cumulative weight loss on a single sample. If the 
single-test-point-per-sample method used with the other materials had 
been used on the U-700 series, over 100 hours of testing time would have 
been necessary to gather sufficient data for a well defined weight loss 
curve. 

Sample 1, aged U-700, was tested for a total of 1080 minutes, as 
shown in Table 5. The sample was removed from the apparatus, washed, 
dried, and weighed at the time intervals shown. The weight column lists 
the sample weight at the end of any test interval; the time column is 
cumulative. The sample was not weighed at the initial damage time. 

Samples 2 and 3 of the solutioned U-700 material were tested to 
initial damage time only. Sample 4, solutioned U-700, was tested for a 
total of 960 minutes, as shown in Table 6. 

REFERENCE 

1. Ripken, J. F., "A Test Rig for Studying Impingement and Cavitation 
Damage," Erosion by Cavitation or Impingement . ASTM Technical 
Publication No. 408, 1967. 
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TABLE 1 


ZINC 


Sample 

V 

(ft/sec) 

Drop 

Diam. 

(mm) 

Impact Exposure 

T. 

(sec) 

Impact 

Rate 

drop/sec 

Initial 

Wt. 

(gms) 

Final 

Wt. 

(gms) 

Comments 

A 1 

1000 

.866 

195.1 

161.5 

.9379 

.9222 

2 

Crater Ridge 

B 

1000 

.866 

135.8 

161.5 

.9490 

.9380 


C 

1000 

.866 

75.7 

161.5 

.9641 

.9600 

Crater Ridge 

2 

1000 

.866 

15.2 

161.5 

.9663 

.9660 

Crater Ridge 

4 

1000 

.866 

30.1 

161.5 

.9724 

.9713 

Crater Ridge 

5 

1000 

.866 

60.3 

161.5 

1.0144 

1.0128 

Crater Ridge 

7 

1000 

.866 

120.4 

161.5 

.9844 

.9801 


lffl 

1000 

.866 

240.2 

161.5 

1.0148 

.9899 

Crater Ridge 

12 

1000 

.866 

90.0 

161.5 

.9816 

.9775 

Crater Ridge 

13 

1000 

.866 

152.5 

161.5 

.9731 

.9621 

Drop Drift 

14 

1000 

.866 

182.7 

161.5 

.9696 

.9514 

No. 14 Hit Off 
Center Hole 
Through Target 

15 

1000 

.866 

180.0 

161.5 

.9938 

.9814 

No. 15 Two Holes 
One 180 Sec. 

One 11.9 Sec. 

16 

1000 

.866 

210.5 

161.5 

.9701 

.9512 


17 

1000 

.866 

225.1 

161.5 

.9603 

.9391 


18 

1000 

.866 

257.1 

161.5 

. 9620 

.9388 

Drill Through 

1 

Specimens A, B, 

And C 

were thicker than 

the other 

specimens . 



A ridge of metal formed around the crater; eroded particles appeared to 
detach from this ridge. 
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IRON 


Sample 

V 

(ft/sec) 

Drop 

Diam. 

(mm) 

Impact Exposure 

T. 

(sec) 

Impact 
Rate 
drop/ sec 

Initial 

Wt. 

(gms) 

Final 

Wt. 

(gms) 

Initial 

Damage 

T. 

(sec) 

1 

1000 

.866 

240.1 

161.5 

1.0645 

1.0631 

170 

2 

1000 

.866 

480.1 

161.5 

1.0664 

1.0640 

210 

3 

1000 

.866 

960.1 

161.5 

1.0376 

1.0321 

180 

4 

1000 

.866 

1920.1 

161.5 

1.0625 

1.0533 

190 

5 

1000 

.866 

3840.1 

161.5 

1.0424 

1.0279 

170 

6 

1000 

.866 

1440.0 

161.5 

1.0515 

1 .0440 

185 

7 

1000 

.866 

2880.2 

161.5 

1.0738 

1.0614 

297 

8 

1000 

.866 

6000.1 

161.5 

1.0634 

1.0414 

197 

10 

1000 

.866 

600.2 

161.5 

1.0604 

1.0556 

121 

11 

1000 

.866 

802.1 

161.5 

1.0482 

1.0425 

157 

12 

1000 

.866 

1207.2 

161.5 

1.0566 

1.0509 

200 

13 

1000 

.866 

5000.0 

161.5 

1.0362 

1.0154 

245 

15 

1000 

.866 

701.1 

161.5 

1.0432 

1.0377 

75 


Comments 


Crater Ridge 


1 


A ridge of metal formed around the crater 



TABLE 3 


NICKEL 


Sample 

V 

(ft/sec) 

Drop 

Diam. 

(mm) 

Impact Exposure 

T. 

(sec) 

Impact 

Rate 

drop/ sec 

Initial 

Wt. 

(gins) 

Final 

Wt. 

(gms) 

Initial 

Damage 

(sec) 

Comments 

1 

1000 

.866 

14.6 

161.5 

1.2242 

1.2233 

10 

Crater Ridge^ 

2 

1000 

.866 

30.4 

161.5 

1.2050 

1.2048 

25 


3 

1000 

.866 

90.1 

161.5 

1.2087 

1.2074 

40 


4 

1000 

.866 

36.0 

161.5 

1.2172 

1.2159 

15 


5 

1000 

.866 

61.9 

161.5 

1.2135 

1.2111 



6 

1000 

.866 

150.2 

161.5 

1.2400 

1.2369 



7 

1000 

.866 

302.3 

161.5 

1.2404 

1.2341 



8 

1000 

.866 

600.2 

161.5 

1.2533 

1.2419 


Crater Ridge 

9 

1000 

.866 

900.0 

161.5 

1.2457 

1.2309 



10 

1000 

.866 

1500.0 

161.5 

1.1957 

1.1760 



11 

1000 

.866 

1200.2 

161.5 

1.2530 

1.2391 



12 

1000 

.866 

1800.0 

161.5 

1.2448 

1.2275 


Crater Ridge 

13 

1000 

.866 

2100.2 

161.5 

1.2524 

1.2318 



14 

1000 

.866 

2400,0 

161.5 

1.2413 

1.2207 



15 

1000 

.866 

2700.7 

161.5 

1.2064 

1.1792 


Footnote 2 


A ridge of metal formed around the crater. 

The eroded hole was oblong which suggests that there was some drift in the point 
of impact of the drop. 
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TABLE 4 


TANTALUM 


Sample 

V 

(ft/ sec) 

Drop 

Diam. 

(mm) 

Impact Exposure 

T. 

(sec) 

Impact 

Rate 

drop/sec 

Initial 

Wt. 

(gms) 

Final 

Wt. 

(gms) 

Initial 

Damage 

(sec) 

Comments 

1 

1000 

.866 

348.7 

161.5 

2.3510 

2.3510 

320 

Crater Ridge 

2 

1000 

.866 

6900.0 

161.5 

2.2779 

2.2286 

270 

Drop Scatter 
3 Drop Diam. 

3 

1000 

.866 

1206.3 

161.5 

2.3154 

2.3055 

230 

Crater Ridge 

4 

1000 

.866 

6000.2 

161.5 

2.2502 

2.2186 

250 


5 

1000 

.866 

2419.6 

161.5 

2.2475 

2. 2267 

170 


6 

1000 

.866 

4799.2 

161.5 

2.2710 

2.2423 

230 


8 

1000 

.866 

600.3 

161.5 

2.3456 

2.3410 

277 

Crater Ridge 

9 

1000 

.866 

9600.0 

161.5 

2.2915 

2.2493 

230 

Hole Thru Target 

11 

1000 

.866 

3250.1 

161.5 

2.2068 

2.1810 

250 


12 

1000 

.866 

2000.2 

161.5 

2.1901 

2.1762 

320 


13 

1000 

.866 

7000.6 

161.5 

2.2161 

2.1805 

- 


14 

1000 

.866 

1500.0 

161 .5 

2.2761 

2.2684 

400 


15 

1000 

.866 

800.2 

161.5 

2.2090 

2.2036 

200 
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TABLE 5 


Sample 
No. 1 


AGED UDIMET 700 


V 

(ft/sec) 

Drop 

Diam 

(mm) 

1000 

.866 

1000 

.866 

1000 

.866 

1000 

.866 

1000 

.866 

1000 

.866 

1000 

.866 

1000 

.866 

1000 

.866 

1000 

.866 

1000 

.866 

1000 

.866 

1000 

.866 

1000 

.866 


Impact Exposure 
T. 

(sec) 

00000 

8100.0 

9000.0 

12616.5 

14416.5 

18016.5 

21616.5 

25216.5 

28816.5 

32416.6 
35999.9 

43200.0 

50400.0 
57600.2 

64800.0 


Impact 
Rate 
drop/ sec 

Weight 

(gms) 

- 

1.1004 

161.5 


161.5 

1.1001 

161.5 

1.0956 

161.5 

1.0927 

161.5 

1.0898 

161.5 

1.0874 

161.5 

1.0854 

161.5 

1.0834 

161.5 

1.0820 

161.5 

1.0806 

161.5 

1.0775 

161.5 

1.0761 

161.5 

1.0742 

161.5 

1.0724 


Comments 


Initial Damage 
(Approximate) 


Hole Thru Target 
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TABLE 6 


SOLUTIONED UDIMET 700 


Sample 

V 

(ft/sec) 

Drop 

Diam. 

(mm) 

Impact Exposure 
T. 

(sec) 

Initial 

Wt. 

(gms) 

Final 

Wt. 

(gms) 

2 

1000 

.866 

3661.2 

1.1626 

1.1623 

3 

1000 

.866 

7200.0 

1.1515 

1.1513 


Sample 
No. 4 

V 

(ft/sec) 

Drop 

Diam. 

(mm) 

Impact Exposure 

T. 

(sec) 

Impact 

Rate 

drop/sec 

Weight 

(gms) 




0000.0 

- 

1.0731 


1000 

.866 

7265.0 

161.5 



1000 

.866 

9065.8 

161.5 

1.0707 


1000 

.866 

10866.2 

161.5 

1.0669 


1000 

.866 

12666.2 

161.5 

1.0658 


1000 

.866 

16266.2 

161.5 

1.0632 


1000 

.866 

19866.3 

161.5 

1.0617 


1000 

.866 

27066.4 

161.5 

1.0586 


1000 

.866 

33553.2 

161.5 

1.0544 


1000 

.866 

41405.7 

161 .5 

1.0500 


1000 

.866 

49605.7 

161.5 

1.0466 


1000 

.866 

57602.9 

161.5 

1.0454 


Comments 


Comments 


Initial Damage 


Hole Thru Target 
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APPENDIX C 


DESCRIPTION OF THE REQUIRED COMPUTER PROGRAM 


by Dr. J. K. Casey, Applied Mathematician, 
Computations Operation, General Electric Company, Evendale, Ohio 


The required program is for the function 


n 

T (n) = ^ v(x)P (x,n) 
x = o 

where 

„ ^ n n! x ^n-x 

P <X ’ n) = “x')— ” P (1_P) 

and 


v (x) = 0 
= R X 

= R]^ + R2 
= R^ + R2 + R3 
= R x + R 2 + R 3 + 


for 0 
for N c 
for N 0 + N x 
for N 0 + N x + N 2 
R 4 for N 0 + N x + N 2 + N3 


< x < N d 

< x < N Q + N-^ 

< x < N 0 + Nx + N 2 

< x < N 0 + Ni + N 2 + N 3 

< x < N d + Njl + N 2 + N 3 + N4 


etc. 

The function v(x) is a step function; the quantities are the 
length of the treads and the R^ are heights of the risers. 

The graph of v(x) should, therefore, have the appearance indicated 
below. v ( x ) 



N 0 N 0 N o +N 1 N 0 +N! 

+N^ +N2 +N 3 +N 3 
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Tread lengths and riser heights beyond those specified in input are 
assumed to equal the last values that were input. Thus if N Q , ... N4 
and R^, ... R4 were input, then N5, Ng, ... are taken equal to N4 and 
R5, Rg, ... are taken equal to R4. Provision should be made for at 
least 7 tread lengths and 6 riser heights to be supplied in the input. 

B. COMPUTER PROGRAM 

by M. A. Cummings, Manager, Mathematical Applications 
Computations Operation, General Electric Company, Evendale, Ohio 

A computer program for the function T(n) was written in double 
precision Fortran IV. The program was based on an earlier version in 
the basic language by Lars H. Sjodahl. A printout of the program is 
given below. 

The input data for the program are: 
probability, p 

counting rule, N 0 , N^, N2, N3, N4, etc. 
volumes, — , R^, R2* ^3» &4» etc. 

number of impacts, n = to n = Step n = 

If the program is run with R^ = R2 = R3 = etc. = 1 , the printout is 
in layers of fragments removed at the specified values of impacts, n. 

If the program is run with the specified values of the R^, the printout 
is accumulated volume loss at the specified values of impacts, n. 

Runs with use of this program can be obtained by applying to 
Mr. M. A. Cummings, Mail Stop H8, General Electric Company, Evendale, 
Ohio. 



DIMENSION NN ( 201 ) *RR ( 200 ) 

DOUBLE PRECISION PO . N , PI ,D »L , X , R 1 , R2 , T ♦ S * P . XO * X 1 , P *D I * R * 
2 NDELTA.NMAX .NN.RR .30UND1 »B0UND2 
DOUBLE PRECISION XX 
NAMEL I ST/ I N/PO.N .NN.RR .NDELTA.NMAX 
WRITE (6,10) 

10 FORMAT (1H1.5X.2HP0.10X.1HN.20X.9HFRAGMENTS// ) 

5 CALL FLGE0F(5,IE0F) 

READ (5, IN) 

IF ( I EOF .NE. 0) GO TO 990 
P1=P0/ t l.-PO) 

100 T=0. 

S = 0. 

I NUM= 1 

SB0UND=NN ( 1 ) 

PSTAR=0. 

D = -1 . 

L=N*PO 

IF ( NN ( 1 ) .GE. L) L=NN ( 1 ) 

X 1 =RR ( INUM) 

S0UND1 =NN ( INUM) 

B0UND2 =NN ( INUM)+NN( INU’1'1) 

IL=IFIX(L) 

X= I L 

R1=X* DLOG(PO) 

R2=(N-X)* DLOG(l.-PO) 

R=0. 

DO 200 1 = 1. IL 
XI =FLOAT ( I ) 

200 R=R+ DLOG ( ( N l.-XIl/XIl 

R= DEXP ( R+ ( N-X ) *DLOG ( 1 .-PO ) +X*DLOG ( PO ) ) 

P=R 

420 IF (X.LT.BOUNDl.AND.BOUNDl.GT.NN(l) ) GO TO 440 
XX=DABS(B0UND1-NN(1) ) 

IF (XX .LT. (.000001) .AND. X .LT. B0UND1 ) GO TO 435 
430 IF ( X.GE. B0UND1 .AND. X.LT. B0UND2) GO TO 445 
IF (INUM .NE. 200) GO TO 433 
WRITE (6.432) 

432 FORMAT (17H PARTIAL SOLUTION) 

XX=X1*P 

WRITE (6.434) PO.N.S.XX 

434 FORMAT ( F 12 . 9 . F 1 1 . 0 , 2 F 1 6 . 9 ) 

GO TO 660 

433 I NUM= I NUM+1 
B0UND1 =B0UND2 
B0UND2=B0UND2 ' NN ( I NUM+1 ) 

X1=X1+RR( INUM) 

GO TO 430 

435 D=l. 

GO TO 600 
440 B0UND2=B0UND1 
I NUM= I NUM-1 

BOUND 1=B0UND1-NN ( I NUM+1 ) 

X1=X1-RR( INUM- 1 1 ) 

445 S=S+P*Xl 

PSTAR=PSTAR+P 

480 IF (D .LT. (0.)) GO TO 500 
p=p*(N— X)/(X+1.)*P1 
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GO TO 510 

500 P=P*X/ ( N-X+l • ) /PI 
510 X =X+D 

IF (<P*X1) .GE. (.000001)) GO TO 420 
590 IF (D .EQ. (-1.0)) SBOUND=X 
IF (D.EQ.1.0) U80UND=X 
D= D+2 • 0 

600 IF (S.EQ.0.0) GO TO 610 
T = T + 5 
S = 0. 

610 IF (D.NE.1.0) GO TO 620 
I X= I F I X ( L+l • ) 

X= FLOA T ( IX) 

IL=IFIX(L) 

P = R* ( N-FLOAT (IL) )/ (FLOAT ( IL )+l. )*P0/ ( 1.-P0) 
GO TO 430 

620 WRITE (6.650) PO »N . S80UND . U80UND . T 
650 FORMAT ( F 1 2 . 9 . F 1 1 . 0 . F6 . 0 , F 5 . 0 . F 1 6 . 9 ) 

660 IF (N .GE. NMAX) GO TO 5 
N = N+NDELT A 
GO TO 100 
990 STOP 
END 


NASA-Langley, 1971 17 JT _ 6532 
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